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ABSTRACT 


This  report  describes  a  SAR  (Synthetic  Aperture  Radar)  image 
quality  analysis  and  an  image  quality  model  incorporating  the 
conclusions  of  the  image  quality  analysis  into  a  SAR  system  digi¬ 
tal  simulation.  Hie  objective  of  the  model  and  its  development 
program  is  to  provide  a  fundamental  approach  for  SAR  image 
quality  effects  for  specific  SAR  systems  (parameter  combinations 
performing  specific  missions).  The  model  thus  provides  a 
relatively  economical  prediction  of  the  information  which  can 
be  extracted  from  a  specified  radar-mission  configuration. 

The  model  output  is  a  collection  of  image  quality  assess¬ 
ment  functions,  in  particular  the  TQF  (Threshold  Quality 
Factor).  Other  parameters,  which  indicate  the  performance 
of  the  SAR  system,  are  available  from  the  digital  simulation. 

These  include  the  point  target  response  (both  rar>ge  and  azimuth) 
along  with  their  locations  for  a  preselected  target  field 
configuration. 


FOREWORD 


This  study  was  :oi.uucted  by  the  Environmental  Research  Institute 
of  Michigan  beginning  in  July  1972.  The  work,  under  contract  F33615- 
72-C-1996,  was  sponsored  by  the  Air  Force  Avionics  Laboratory, 

Air  Force  Systems  Command,  at  Wright-Patterson  Air  Force  Base, 
Dayton,  Ohio.  The  project  engineer  was  Sgt.  T.  F.  Stipulkosky. 

The  project  was  directed  at  ERIM  by  Dr.  L.  J.  Porcello  of  the 
Institute's  Radar  and  Optics  Division.  The  principal  investigators 
were  Dr.  W.  D.  Hall  and  Dr.  R.  H,  Mitchel.  Other  contributors 
to  the  project  included  Mr.  M.  Evans  and  Mr.  B.  Friberg  who  are 
responsible  for  the  AIRMOD  module,  Mr.  H.  Doss  who  generated 
the  RSMOD,  and  Dr.  D.  Ausherman  who  devised  the  IRMOD  module. 
The  image  quality  analysis  and  IQMOD  development  were  performed 
by  Dr.  R.  Mitchel. 

This  final  report  consists  of  two  volumes.  Volume  I  is  CONFI¬ 
DENTIAL  and  describes  the  image  quality  analysis  and  its  conclusions 
along  with  an  overview  of  the  SAR  simulation.  The  second  volume  is 
Unclassified  and  contains  detailed  descriptions  of  each  of  the  modules 
composing  the  simulation. 

Details  of  the  technical  effort  for  the  period  1  July  -  31  December 
1972  are  contained  in  the  interim  report,  "Image  Quality  Analysis,  " 
ERIM  Report  193900-1-P,  Ralph  H.  Mitchel,  July  1973. 


TABLE  OF  CONTENTS 


Introduction .  1 

5.0.  AIRMOD  Module  Description .  3 

5. 1.  Mode,  description . . .  3 

5.2.  Computer  Program . , . .  12 

6.0.  Radar,  Recorder  and  Processor  Simulation .  18 

6.1.  Introduction  and  Background .  18 

6.2.  Main  Program  and  Input  Models .  24 

6.2.1.  Targets . 24 

6.2.2.  Aircraft  Parameters .  28 

6.2.3.  Transmitter-Receiver .  29 

6.2.4.  Motion  Compensation .  34 

6.2.6.  Processor .  42 

6.2.7.  Displays .  45 

6.3.  Signal  Generation .  46 

6.3.1.  Timer .  47 

6.3.2.  Table .  51 

6.3.3.  Point .  51 

6.3.4.  MOCOMP .  54 

6.3.5.  AMPR .  54 

6.3.6.  CHIRPX .  60 


v 


6.3.7.  LIMIT .  60 

6.3.8.  MTF's .  63 

6.4.  Processor . . . .  64 

6.4.1.  Focusing  for  Signal  Histories .  64 

6.4.2.  Optimization  Choices. . . .  73 

6.4.3.  Image  Location . . .  76 

6.4.4.  Geometrical  Distortion .  77 

6.4.5.  Processor  Subroutines .  79 

6.5.  Output .  80 

6.6.  Conclusions .  95 

7.0.  Impulse  Response  Generation .  97 

7.1.  Introduction .  97 

7.2.  Theory  of  Operation .  .  . . .  97 

7.2.1.  Received  Signal  Generation .  100 

7.2.2.  Reference  Function  Generation .  104 

7.2.3.  Performing  the  Convolution. .  106 

7.3.  Target  Response  Software  Description . 109 

7.4.  Processing  Examples . 119 

7.4.1.  Demonstration  of  Motion  Errors .  132 

7.4.2.  Demonstration  of  Amplitude  Weighting  Effects  .  .  .  141 

7.5.  Summary . 156 


Appendix 


157 


LIST  OF  ILLUSTRATIONS 


1.  Single  Channel  of  Inertial  Platform .  13 

2.  Vertical  Channel .  14 

3.  Beam  Velocity  Computer .  15 

4.  Ground  Speed  Channel .  16 

5.  Model  Block  Diagram .  19 

6.  Radar  Beam  and  Target  Geometry . 21 

7.  Main  Program  Generalized  Flow  Chart . 25 

8.  Amplitude  of  Target  Return  vs.  Time .  32 

9.  Forward  Squint  Angle . 37 

10.  SIGGEN  Generalized  Flow  Diagram . 48 

11.  Antenna  Beam  Width . 50 

12.  Subroutine  Table . 52 

13.  Antenna  Gain  -  Vertical  Position . 57 

2 

14.  Horizontal  Beam  Structure  (sine) .  58 

15.  Amplifier-Limiter  Function . 62 

16.  Range  Response  MTF . 65 

17.  Azimuth  Response  MTF . 66 

18.  Modulation  Transfer  Function . 67 

19.  Recorder  Limiting  Function . 68 


UST  OF  ILLUSTRATIONS 
(continued) 


20.  Sample  Plots  FS  vs.  Z  .  .  7  i 

21.  Corrections  for  Squint  and  Distortion .  78 

22.  PROCES  Generalized  Flow  Diagram .  81 

23.  Signal  Pulse  Amplitude  as  a  Function  of  Film  Position.  .  .  .  121 

24.  Doppler  Frequency  of  Return  as  a  Function  of  Azimuth 

Position .  123 

25.  Signal  Phase  as  a  Function  of  Azimuth  Position .  124 

26.  One-Dimensional  Zone  Plate  Section  Constructed  Using  Signal 

Phase  . 125 

27.  Quadratic  Reference  Phase  as  a  Function  of  Azimuth  Position  127 

28.  Point  Target  Response . 128 

29.  An  Example  of  the  Three  Plot  Options  Associated  with  the 

Point  Target  Response .  129 

30.  Idealized  System  Impulse  Response .  131 

31.  Far  Range,  No  Motion  Point  Target  Response .  134 

32.  Mid  Range,  No  Motion  Point  Target  Response .  135 

33.  Near  Range,  No  Motion  Point  Target  Response .  136 

34.  Far  Range  Point  Target  Response  with  Motion  Error .  ...  137 

35.  Mid  Range  Point  Target  Response  with  Motion  Error.  .  .  .  138 

36.  Near  Range  Point  Target  Response  with  Motion  Error  ...  139 

37.  Uniform  Amplitude  Weighting  for  Returned  Signal  Pulse  .  .  142 

38.  Sin  x/x  Amplitude  Weighting  (out  to  -3  dB)  for  Returned 

Signal  Pulse .  143 


viii 


LIST  OF  ILLUSTRATIONS 
(continued) 


39.  Uniform  Amplitude  Weighting  Far  Range  Response .  145 

40.  Uniform  Amplitude  Weighting  Mid  Range  Response .  146 

41.  Uniform  Amplitude  Weighting  Near  Range  Response .  147 

42.  Non-Uniform  Weighting  Far  Range  Response .  148 

43.  Non-Uniform  Weighting  Mid  Range  Response .  149 

44.  Non-Uniform  Weighting  Near  Range  Response .  150 

45.  Three  Sidelobe  Si.i  x/x  Antenna  Weighted  Doppler  History  .  .  .  152 

46.  Far  Range  Target  Response .  153 

47.  Mid  Range  Target  Response . . .  154 

48.  Near  Range  Target  Response .  155 


{ 

t 

} 


INTRODUCTION 

This  is  the  second  volume  of  a  report  describing  a  SAR  (Synthetic  Aperture 
Radar)  image  quality  analysis  and  an  SAR  image  quality  model  incorporating 
the  conclusions  of  the  image  quality  analysis  into  a  SAR  digital  simulation. 

The  objective  of  the  model  and  its  development  program  is  to  provide  a 
fundamental  approach  for  SAR  image  quality  analysis  via  a  model  containing 
sufficient  sophistication  to  quantify  image  quality  effects  for  specific  SAR 
\  systems  (parameter  combinations  performing  specific  missions) .  The  model 

r  thus  provides  a  relatively  economical  prediction  of  the  information  which  can 

|  be  extracted  from  a  specified  radar-mission  configuration. 

t 

|  The  approach  used  to  define  an  image  quality  figure  of  merit  and  the 

description  of  the  resulting  image  quality  model  is  the  subject  of  this  report. 

The  report  is  written  in  two  volumes.  The  first  volume,  describes  the  image 
quality  analysis  and  its  conclusions  along  with  an  overview  of  the  SAR  simulation. 
The  second  volume  contains  detailed  descriptions  of  each  of  the  modules  composing 
the  simulation. 

Sections  5,  6,  and  7  in  this  volume  describe  the  details  of  the  simulation 
modules.  Section  5  discusses  the  AIRMOD  module,  used  to  compute  the  airframe 
-  dynamics.  Section  6  describes  the  RSMOD  module  which  simulates  the  radar, 

recorder  (storage),  and  p ’ocessor.  Section  7  details  the  theory  and  implementation 
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of  the  IRMOD  module  which  computes  the  impulse  response  characteristics 
of  the  entire  SAR  system. 


5.0  AIRMOD  MODULE  DESCRIPTION 

The  aircraft  model  (laceral  channel  only)  was 
originally  programmed  on  the  IBM  360  digital  computer 
to  simulate  the  University  of  Michigan  C-131  aircraft. 
Later  simulations  included  the  autopilot,  a  simplified 
doppler  navigator,  and  a  single  axis  inertial  platform. 

The  current  model  is  an  expansion  of  the  previous 
program  that  includes  the  longitudinal  and  vertical 
dynamics,  a  r.iore  complex  autopilot  and  a  turn  command 
generator.  The  aircraft  can  be  commanded  to  fly  a 
straight  line  or  a  complex  three  dimensional  flight 
path.  The  model  for  the  inertial  platform  includes  the 
roll,  pitch  and  yaw  axes,  the  dynamics  for  three  accelero 
meters  mounted  on  the  X,  Y,  Z  axes,  and  platform  tilt 
errors.  The  model  has  the  built  in  flexibility  to  simu¬ 
late  most  modern  aircraft  plus  the  flexibility  to  simu¬ 
late  many  aircraft  radar  configurations. 

5.1  MODEL  DESCRIPTION 

The  aircraft  is  considered  a  rigid  body,  and  the 
path  of  the  aircraft  is  determined  by  the  aerodynamic 
forces,  the  angular  acceleration  about  the  three  ortho¬ 
gonal  axis,  the  inertia  characteristic,  the  propulsive 
force,  and  the  air  density.  The  eauations  of  motion 


are  derived  from  the  basic  force  and  moment  equations 
for  each  degree  of  freedom.  These  force  equations  are 
then  written  in  partial  differential  form.  Closed 
loop  operation  of  the  aircraft  is  maintained  through  the 
control  of  the  aileron,  rudder  and  elevator.  Control 
surface  deflections  produce  the  desirable  moments  along 
the  principal  axes.  These  deflections  also  generate 
cross  moments  and  forces  along  the  principal  axes  that 
do  not  always  contribute  positively  to  the  control  of 
the  aircraft. 

The  model  simulates  the  aircraft  using  nine  degrees 
of  freedom,  the  lateral,  vertical,  and  along  track 
accelerations,  the  roll,  pitch,  and  ya»  rotations,  and 
the  three  control- surface  deflections.  The  six  dif¬ 
ferential  equations  can  be  divided  into  two  groups  of 
three.  The  first  three  equations  are  the  lateral  side 
force  equation,  the  roll  equation,  and  the  yaw  equation. 
These  equations  are  written  in  terms  of  the  variables 
6  (drift  angle) ,  ^  (yaw  angle) ,  4>  (roll  angle) ,  6a 
(aileron  deflection  angle),  and  6r  (rudder  deflection 
angle) .  The  second  three  equations  are  the  longitudinal 
force  equation,  the  vertical  force  equation,  and  the 
pitch  angle  equation.  The  variables  are  V  (aircraft 
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velocity) ,  a  (angle  of  attack)  ,  9  (pitdh  angle)  ,  and 
6e  (elevator  angle) .  The  nine  degree  of  freedom 
equations  are  written  below: 

1)  a  *  *  -  If  sin  *  ♦  \  £™  ,  c  „ 

*  W') 

2)  ; ,  iq^a.  [ctB#  *  ctte«»  *  cMr*r 

*  TT  (Clp*  *  c*r*)] 

S)  ; .  .  Wa  ,  Wr 

*  "T1  (cnp*  *  cnr*)] 

l)  •  *  [Cma  4“  *  We  +  W  (CmiS  *  Cme5)] 


C£B3  +  CUa6a  +  CUr6r 


C  „B  +  C  »  6a  +  C  ,  6r 
nB  n6a  n6r 


•6  .  SIMJ1  |cma  ,a  ♦  Cm6e6e  ♦  ^ 


5)  mV (a  - 


a  -  9)  »  lpV2S  T — V-  + 

[kv 


CLa  Aa  +  CL6eA5e  “M 


7  Py2s 


T~ZT:  +  CLa  Aa  +  CL6eA6e 


*VoS 


r  (2.5A)  -  W  (9  -  a) 
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The  static  stability  derivatives  C^g ,  C^g,  Cng  etc., 
and  the  dynamic  stability  derivatives  ,  C^,  and 
Cnr  etc.  are  some  of  the  parameters  that  describe  the 
aircraft  under  investigation.  These  plus  other  constants 
such  as  the  mass,  the  moment  of  inertia  about  each 
axis,  etc.  are  calculated,  measured,  and/or  determined 
from  wind  tunnel  tests  and  then  substituted  in  the 
above  equations.  Host  conventional  aircraft  can  then 
be  readily  simulated  using  the  above  equations.  The 
equ'J^ions  can  then  be  solved  by  using  classical  techniques, 
a  rather^ormidable  task,  or  the  equations  can  be  pro 
grammed  intfLa  computer.  There  are  several  second  order 
and  nonlinear 'effects  that  have  been  assumed  negligible. 

The  control  surface  dynamics  are  approximated  by  second 
order  differential  equations,  and  the  resulting  errors 
are  assumed  to  have  little  influence  on  the  radar,  but 


may  affect  the  aircraft  handling  qualities.  A  selection 

from  any  of  the  above  variables  can  be  listed  as  a 

function  of  time  in  the  output  file. 

The  aircraft  motions  are  available  from  the  system 

dynamics  and  are  related  to  a  system  of  coordinates.  A 

flight  path  coordinate  system  is  labeled  -  along  track, 

-  lateral  direction,  ZQ  -  vertical  direction,  and  the 

fixed  axes  are  labeled  X  ,  Y  ,  Z  .  The  transformation 

ooo 

matrix  between  the  two  can  then  be  written. 


and  the  inverse  matrix  written 


An  aircraft  oriented  reference  frame  is  also 

written  labeled  X  ,  Y  ,  Z  .  and  the  transformation  mat- 

a  a  a 

rix  is  written  with  respect  to  the  "f"  reference  system. 


The  matrix  is  also  written  for  the  transformation  between 
the  "a"  and  the  "o"  axis  by 

[A,,om]  =  [Tof]  [A"f"] 

The  velocity  matrix  for  the  "o"  reference  frame  is 
written  by  integrating  the  acceleration  signals 

[v,,o"] +  [/]  M 

and  a  displacement  matrix  is  written  by  a  second 
integration 

M =  [/]  [v,,o’] 


A  turn  command  generator  is  integrated  in  the 
model  that  has  the  capability  of  introducing  simple 
or  complex  turns  in  the  lateral  axis  and/or  the  ver¬ 
tical  axis. 

The  model  is  commanded  to  a  straight  flight  path 
(subject  to  atmospheric  turbulence)  when  the  turn 
generator  is  off.  The  turn  generator  is  switched  "ON" 
by  selecting  a  non  zero  value  for  one  or  more  of  the 
constants  in  the  following  expression 

a  =  +  I^t  +  K^t^  +  sin  (W^t  +  <f>^)  +  sin  (l^t  + 

As  an  example,  a  selection  of  a  value  initiates  a 
constant  acceleration  command  at  time  t  =  0.  The  air¬ 
craft  control  surface  is  deflected  and  a  transient  is 
initiated  in  the  aircraft.  A  settling  period  is  required 
before  the  acceleration  error  signal  reaches  a  steady 
state  value.  The  model  also  has  the  capability  of 
starting  at  t  -  0  in  a  settled  condition,  i.e. ,  the 
autopilot  error  signal  maintains  a  steady  state  value. 

The  complex  turn  command  also  includes  a  second  expression 
that  is  used  to  control  the  vertical  channel,  and  is 
independent  from  the  lateral  channel. 

The  doppler  navigator  is  a  radar  system  that  measures 
aircraft  ground-speed  and  drift  angle.  The  model 
approximates  the  doppler  by  a  first  order  system  with 
error  inputs  (El)  added.  A  perfect  doppler  navigator 
is  also  modeled  to  show  the  possible  improvement  in 
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the  overall  system  performance.  In  the  "base"  system 
the  doppler  drift  angle  is  combined  with  azimuth  angle 
information,  and  this  signal  is  sent  to  a  steering  gyro. 
The  sensitive  axis  is  mounted  to  detect  yaw  motion. 

(The  terminology  "base"  system  refers  to  an  aircraft  - 
inertial  navigation  subsystem  that  has  been  used  in  a 
previous  model.  A  comparison  can  then  be  made  between 
the  new  and  old  models  to  confirm  programming  accuracy.) 
The  gyro  output  is  then  sent  to  the  autopilot  rudder 
channel. 


0  -  drift  angle 
t|;  -  azimuth  angle 
K  -  arbitrary  gain 

The  "base"  autopilot  is  very  simply  approximated 
by  a  displacement  channel  and  a  rate  channel  for  each 
control  surface.  The  servo  errors,  assumed  to  be 
small,  were  not  included  for  the  base  study.  The  "base' 
autopilot  is  incorporated  in  the  autopilot  diagrams 
shown  in  the  figures.  The  autopilot  block  diagrams 


for  the  three  channels  are  given  in  Figures  2-4  in 
the  computer  program  section  5.2.  The  aileron 
channel  has  three  sets  of  inputs,  the  true  values 
(acceleration,  velocity,  displacement,  roll  and  roll 
•  rate),  the  inertially  computed  values  (acceleration, 

velocity,  and  displacement).  A  vertical  gyro  is  included 
that  is  oriented  to  the  local  vertical. 


— 

Y 

-* 

1 

a 

A. 
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The  gyro  pendulum,  being  sensitive  to  lateral  accelera¬ 
tions,  is  torqued  off  level  during  a  turn,  and  a  settling 
period  is  required  to  restore  the  aircraft  to  a  straight 
and  level  flight  path.  The  autopilot  rudder  channel  also 
includes  three  sets  of  inputs,  the  ideal  signals  (accelera¬ 
tion,  velocity,  displacement,  and  azimuth  rate),  the  com¬ 
ponent  generated  signals  (acceleration,  velocity,  dis¬ 
placement,  and  azimuth),  and  the  command  signals  (accel¬ 
eration,  velocity  and  displacement).  The  autopilot  ele¬ 
vator  channel  also  includes  the  three  sets  of  inputs. 

In  addition,  a  pitch  gyro,  a  ptich  rate  gyro,  a  barometric 
altimeter,  and  a  barometric  rate  unit  are  approximated 
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in  the  model.  The  above  signals  are  combined  according 
to  the  gain  settings,  and  sent  to  the  three  servo  sys¬ 
tems  to  control  the  elevator,  rudder  and  aileron. 

The  inertial-platform  (INS)  model  is  represented 
by  four  subsystems:  the  lateral  channel,  the  longi¬ 
tudinal  channel,  the  azimuth  channel,  and  the  vertical 
channel.  The  first  two  channels  are  functionally 
identical.  Figure  1  is  a  diagram  for  a  single  channel 
inertial  platform,  that  includes  the  accelerometer 
dynamics  and  the  platform  tilt  errors  that  are  a  fun¬ 
ction  of  angle  inputs.  Accelerations  are  sensed  by 
the  accelerometer  mounted  on  the  inner  gimbal,  and 
integrated.  The  velocity  signals  are  either  subtracted 
from  the  doppler  navigator  outputs  or  used  directly  to 
torque  the  gyro,  depending  upon  the  choice  of  doppler 
damped  or  Schuler  mode  of  operation.  The  velocity 
signals  are  then  integrated  to  produce  displacement 
information.  The  model  also  computes  the  change  in  the 
local  vertical  as  the  aircraft  travels  over  the  surface 
of  the  earth.  Gain  control  is  provided  to  select  any 
degree  of  doppler  damping.  The  platform  azimuth  channel 
senses  the  azimuth  angle  information  by  means  of  a  vertical 
gyro.  The  gyro  torqueing  signal  is  derived  from  the  com¬ 
puted  velocity,  earth  rate,  and  latitude  information. 

The  platform  vertical  acceleration,  the  barometric 
altitude,  and  altitude  rate  are  the  inputs  to  the  vertical 
channel  in  Figure  2.  Barometric  and  barometric  rate 
inputs  are  provided  for  damping  purposes.  The  outputs 
are  sent  to  the  autopilot  elevator  control  channel. 
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Accelerations  derived  from  the  platform  are  sent  to 
the  beam-velocity  computer  (Figure  3) .  In  this  subsystem, 
the  aircraft  velocity  is  computed  along  a  predetermined 
depression  angle  and  squint  angle.  The  feedback  loop  is 
the  sum  of  an  integrator,  a  filter,  and  an  arbitrary  gain, 
and  is  subtracted  from  the  acceleration  input.  The  com¬ 
puter  can  be  operated  :n  either  a  first,  second,  or  third 
order  damped  mode  with  any  choice  of  loop  gain  and  damping. 
The  second  channel  receives  the  true  acceleration  signals, 
resolves  them  along  the  predetermined  angle,  and  inte¬ 
grates  to  obtain  the  true  velocity:  Three  basic  outputs 

are  available,  V,  ,  V  ,  and  V,  . 

’  bm’  t  be 

The  ground  speed  channel  in  Figure  4  computes  the 
ground  speed  using  the  combined  inertial-platform-derived 
ground  track  acceleration  signals  axm,  a^  and  the  doppler- 
derived  ground  speed  V  .  The  combined  output  velocity 

O 

contains  both  doppler  (low  frequency)  and  platform  (high 
frequency)  data.  The  error  computer  channel,  indicated 
by  the  dashed  lines,  accepts  the  true  acceleration  signals, 
resolves,  integrates  and  compares  with  the  computed  velocity 
to  produce  the  velocity  error  signal  V  . 

The  error  computer  is  not  a  separate  unit,  but  is 
integrated  throughout  the  model.  The  error  computer  com¬ 
bines  the  true  values  of  the  system  variables  with  the 
component  derived  signals.  These  signals,  representing  the 
subsystem  errors  as  a  function  of  time,  are  stored  in  a 
file  and  can  be  printed  for  documentation. 
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Velocity  Computer 
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5.2  COMPUTER  PROGRAM 

The  AIRMOD  program  is  written  in  the  CSMP  language 
that  is  compatible  with  the  IBM  360  computer.  The  CSMP 
format  has  many  optional  operating  modes,  that  allows  for 
considerable  flexibility  in  the  programming.  The  dynamic 
section  contains  the  bulk  of  the  program:  the  set  of  dif¬ 
ferential  equations  that  relates  to  the  system  under  investi¬ 
gation.  The  use  of  the  initial  section  is  optional  but  an 
improvement  in  computer  efficiency  is  obtained.  Other 
sections  include  the  initial  conditions,  the  timing 
commands,  and  the  output  specifications. 

The  output  file  contains  the  time  histories  of  the 
many  system  variables  that  are  used  as  inputs  for  the 
MSTRIP  computation.  The  general  organization  of  this  simu¬ 
lation  allows  the  operator  to  select  a  variety  of  configura¬ 
tions,  gains,  and  operational  parameters  suitable  for  exer¬ 
cising  an  aircraft  in  a  multitude  of  SAR  modes.  The 
importance  of  this  sumulation  is  the  ability  to  investigate 
the  results  of  hardware  parameter  selection  (e.g. ,  platform 
sophisitcation)  and  effects  of  aircraft  maneuvers  on  the 
overall  SAR  mapping  performance. 
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6.0  RADAR,  RECORDER  AND  PROCESSOR  SIMULATION 


6.1  Introduction  and  Background 


A  digital  computer  simulation  is  described  in  this  section. 

Its  basis  is  a  mathematical  model  in  which  important  hardware 
components  are  represented  and  associated  together  to  form  a  system 
simulation.  Figure  5  illustrates  hardware  components  and  system 
flow.  Inputs  and  outputs  to  the  simulation  are  flexible  to  provide 
convenience  and  generality.  Input  parameter  values  can  be  chosen 
to  reflect  many  different  system  characteristics  and  then  easily 
varied  from  run  to  run  to  facilitate  parameter  studies  as  required. 
Outputs  include  measures  of  system  performance  [7.0]  and  inter¬ 
mediate  results  that  are  useful  as  inputs  for  additional  analysis 
[2.0].  An  original  form  of  this  simulation  was  named  MSTRIP 
(strip-map  radar  simulation).  The  present  model  has  been  named 
RSMOD  (Radar,  recorder  and  processor  simulation). 

Aircraft  motion  and  target  motion  can  be  either  inputted  as 
tables  from  external  sources  such  as  AIRMOD,  discussed  in  a  pre¬ 
vious  section,  and  logs  from  aircraft  flights,  or  generated 
internally  by  function  generators.  Motion  compensation  (e.g., 
phase  snifter,  sampling  rates,  trigger  delays,  clutter  lock) 
may  be  generated  internally  or  can  be  externally  generated  by  a 
flight  simulator  such  as  AIRMOD. 

Radar  parameters  are  set  up  as  inputs  of  wavelength,  nominal 
squint  and  depression  angle,  half-power  beam  width,  scanning 
beam  width,  and  looking  beamwidth.  Antenna  patterns  are  included 
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Figure  5.  Model  Block  Diagram 
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by  tabular  Input  and/or  through  the  use  of  the  sine  function 

ttD 

where  x  ■  y-  sin  8  with  D  antenna  aperture,  X  wave¬ 
length,  and  6  as  the  look  angle  from  the  antenna.  The  antenna 
pattern  is  maintained  vertically  and  horizontally,  and  is  directed 
toward  the  target  field  by  the  instantaneous  pointing  angles  of 
the  radar.  Half-power  beam  widths  (BEAM)  are  used  to  set-up  the 
sine  function  for  half-power  beam  inputs.  Scanning  beamwidths 
(BEAMSC)  relate  the  rectangular  beamwidth  over  which  point  targets 
will  be  interrogated  for  chirp  and  doppler  sampling  (i.e.,  gener¬ 
ation  of  signal  history) .  The  looking  beamwidth  (BEAMLK)  is  used 
to  obtain  background  clutter  by  integrating  clutter  over  the 
sensible  beamwidth  and  to  allow  amplitude  distortion  by  target 
returns  entering  through  side  lobes.  BEAM  <  BEAMSC  £  BEAMLK, 
where  equality  is  not  representative.  However,  one  may  input 
equals  if  desired.  Internally  generated  motion  compensation  may 
be  driven  relative  to  BEAM  or  BEAMSC. 

Geometry  of  aircraft,  target  and  radar  beam  is  considered 
in  a  realistic  manner.  Figure  6  illustrates  the  basic  geometry 
in  which  the  radar  is  looking  out  the  right-hand  side  of  an  air¬ 
craft.  Looking  out  the  left  side  is  obtained  by  using  a  squint 
angle  greater  than  90°.  A  squint  of  about  30°  is  illustrated. 

Three  cooridnate  frames  of  reference  are  presented.  The  a-frame  relates 
aircraft  rotation  and  translation,  the  s-frame  relates  target 
motion,  and  the  (x,  y,  z)  frame  is  the  master  frame.  Time  depen¬ 
dent  transformations  are  maintained  to  relate  all  motion  to  the 
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master  frame.  During  a  simulation  run,  aircraft  velocities, 
target  velocities,  and  antenna  pointing  angles  are  used  to  compute 
instantaneous  beam  velocities  from  which  signal  histories,  tem¬ 
poral  and  spatial,  are  computed.  Motion  compensation,  selec¬ 
table  either  from  a  number  of  possible  functions  or  from  tabular 
inputs,  relates  the  temporal  histories  to  the  recorded  spatial 
histories.  Signal  histories  in  doppler  and  chirp  (video)  are 
made  dependent  on  the  carrier  frequency  and  the  chirp  pulse 
that  is  transmitted.  From  the  received  chirp  and  doppler  a  video 
frequency  and  amplitude  history  for  each  of  n  (n  «  21  is  a  present 
programming  limit  which  can  be  changed  by  recompiling  the  program) 
point  targets.  This  history  reflects  signal  amplitude,  background 
clutter,  and  noise  level  as  functions  of  radar  gain,  RF  gain  and 
limiting,  video  gain  and  limiting,  modulation  transfer  functions 
applied  during  signal  recording,  and  the  transmitted  power  of  the 
radar.  A  target  pattern  of  cross-section  is  provided  for  each 
point  target. 

Output  from  simulation  runs  can  be  varied  also.  Tables  of 
doppler  and  video  history,  and  binary  coded  signal  histories  for  direct 
input  to  impulse  response  calculations  can  be  chosen  for  outputs. 

A  linear  fit  processor  may  be  applied  to  compress  (focus)  video 
and  doppler  histories  into  range  and  azimuth  positions,  and  range, 
azimuth,  and  RMS  (Range  and  Azimuth)  resolutions.  The  ratio  of  input 
distance  to  output  distance,  between  each  pair  of  point  targets,  is 
computed  as  a  measure  of  geometric  distortion.  A  digital  plot  of 
input  images  and  output  images  of  the  point  targets  is  also  available. 


To  conclude  this  Introductory  discussion,  we  observe  that 
AIRMOD  or  real  aircraft  motion  histories  along  with  real  target 
motion  (e.g.,  ships)  can  be  used  as  inputs  to  RSMOD.  That  RSMOD 
is  a  hardware  oriented,  parametric,  representation  if  a  strip- 
map,  synthetic  aperture,  radar.  Mathematical  approximations  are 
used  to  restrict  computation  but  are  chosen  to  provide  a  degree 
of  similitude  which  gives;  simulation  results  that  are  comparable 
with  real  systems.  Although  an  optical  recorder  and  processor  are 
presently  represented,  modular  construction  of  the  computer  mode] 
allow  large  or  small  changes  by  re-programming.  Large  changes 
would  be  represented  by  changing  from  chirp  to  a  different  type 
radar  pulse  and  by  changing  from  one  processor  type  to  another. 

Small  changes  could  involve  any  number  of  function  generators 
such  as  a  chirp  generation,  a  limiter  function,  or  the  sine  func¬ 
tion  generator.  Wherever  possible  function  generators  have  been 
singled  out  as  subroutines  to  allow  easy  transition  to  different 
or  more  complex  forms.  Even  so,  the  present  simulation  is  not 
tied  to  any  particular  hardware.  Through  the  use  of  a  multitude 
of  input  parameters,  many  different  aircraft,  radars,  radar  controls, 
and  other  system  components  are  represented  by  RSMOD  in  its  pre¬ 
sent  form. 

The  remainder  of  this  section  is  devoted  to  a  comprehensive 
discussion  of  simulation  modules. 
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relative  to  a  target  frame  of  reference.  This  target  frame, 
shown  in  Fig.  6  as  the  s-frame,  has  coordinates  relative  to  the 
main  frame  which  may  be  arbitrarily  chosen  in  position  and  in 
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Figure  7.  (a)  Main  Program  Generalized  Flow  Chart 
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angular  orientation.  Each  target  point  may  be  given  an  arbitrary 
maximum  cross-section  with  a  vertical  and  horizontal  beamwidth. 

What  this  means  is  that  each  target  can  reflect  as  if  it  resembled 
a  flat  plate  (very  narrow  beamwidth)  or  as  a  corner  reflector 
(very  wide  beamwidth)  or  anywhere  in  between.  The  cross-section 
inputs  consist  of  a  maximum  cross-section,  the  two  beamwidths,  and 
two  pointing  angles,  one  in  azimuth  and  the  other  in  elevation, 
which  describe  how  the  target  reflects  directionally  relative 
to  the  target  frame.  A  Sine  function  is  then  computed  from  these 
inputs  which  produces  a  cross-section  pattern  with  a  main  lobe 
structure  which  has  vertical  and  horizontal  1/2-power  beam  struc¬ 
ture  corresponding  to  the  input  beamwidths. 

As  noted  earlier,  individual  target  locations  are  input 
relative  to  a  target  frame  of  reference.  Initial  inputs  set 
the  beginning  location  and  orientation  of  the  target  frame. 

Hence,  individual  point  targets  present  a  motion  history  that  is 
directly  related  to  position  in  the  target  frame. 

Target  and  radar  interactions  are  computed  on  a  per  target 
basis.  Amplitude  of  target  returns  are  dependent  on  cross-section 
of  targets  and  the  antenna  pattern  of  the  radar  which  in  turn  are 
aspect  angle  dependent.  Thus,  targets  may  scintillate  to  some 
extent  relative  to  the  interaction  between  two  lobed  patterns. 

High  frequency  scintillation  is  restricted  by  the  computational 
load  inherent  in  a  large  simulation  program.  Of  course,  many  of 
the  overall  effects  of  a  closely  knit  target  complex  can  be  studied 
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through  the  description  of  cross-section  vs.  aspect  of  a  single 
target  center. 

The  limitation  of  having  only  point  targets  with  amplitude, 
angle-dependent,  return  characteristics  represents  a  significant 
simplification  which  makes  the  simulation  tractable  from  a  cost 
and  programming  point  of  view.  Using  previously  acquired  radar 
target  signature  data,  various  target  complexes  and  target  types 
can  be  quite  easily  simulated  by  arranging  the  point  targets  in 
a  suitable  configuration  representative  of  these  selected  target 
types.  This  approach  precludes  the  problem  of  constructing  a 
microwave  scattering  theory  suitably  realistic  to  model  targets 
in  the  complex  world.  This  approach  does  introduce  some  depar¬ 
ture  from  reality  in  that  dispersive  targets  are  not  included, 
but  for  the  problem  we  are  initially  addressing,  it  does  not 
present  insurmountable  departures  from  the  overall  SAR  system 
performance  evaluation  task. 

6.2.2  Aircraft  Parameters 

Aircraft  input  parameters  may  be  used  in  the  AIRMOD  simula¬ 
tion  to  provide  a  time  history  of  aircraft  motion.  Inputs  include: 

a.  initial  airspeed 

b.  altitude 

c.  gross  weight 

d.  stability  coefficients 

e.  control  coefficients 

f.  autopilot  gains 


g.  steering  commands 

h.  wind  velocity 

i.  gust  profiles 

j.  INS  (Inertial  Navigation  System) 

k.  doppler  errors 

l.  trubulence  level 

m.  antenna  stability  errors 

n.  desired  flight  path 

The  AIRMOD  simulation  can  be  used  to  obtain  time-histories 
of  aircraft  altitude,  position,  velocities,  accelerations  or  any 
other  AIRMOD  parameter  and  its  error.  Limitations  to  AIRMOD 
are  the  linear,  rigid-body,  differential  equations  used  for  air¬ 
frame  representations.  Otherwise,  this  aircraft  simulation  pro¬ 
vides  realistic  parameters  with  errors  for  input  to  RSMOD. 

We  note  that  RSMOD  can  be  run  without  AIRMOD  by  limiting 
aircraft  motion  to  prescribed  time-functions  of  antenna  transla¬ 
tion  and  rotations.  Inputs  to  these  functions  are  coefficients 
which  describe  position,  velocity,  magnitude,  frequency  and 
phase  of  a  sinusoidal  term  and  magnitude  of  a  step  function 
term.  Aircraft  motions,  recorded  in  digital  format  from  aircraft 
flight,  may  be  used  as  inputs  also.  In  which  case,  the  only 
difference  from  AIRMOD  inputs  would  be  the  pre-processing  of 
records  into  an  acceptable  table  for  input  to  RSMOD. 

6.2.3  Transmitter-Receiver 

Inputs  to  the  transmitter-receiver  section  are: 
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a.  wavelength  (nominal) 

b.  squint  and  depression  angles  (nominal) 

c.  three  beamwidths  (geometrical) 

d.  antenna  pattern 

e.  chirp  function  with  non-linearities 

f.  transmitter  power 

g.  antenna  gain  (transmitter  and  receiver) 

h.  receiver  front-end  noise 

i.  RF-amplifier  gain  and  limiting 

j .  detector  noise 

k.  video-amplifier  gain  and  limiting 

l.  antenna  stabilization  errors 

m.  PRF 

n.  logical  switches 

o.  background  clutter  levels. 

An  explanation  of  these  radar  parameters  follows.  These 
radar  parameters  are  used  in  the  following  manner.  For  MSTRIP, 
beamwidth  was  specified  as  an  angle  spreading  out  from  the  radar. 
This  pattern  was  aimed  by  a  specification  of  squint  (forward  from 
broadside)  and  depression  (positive  downward  from  the  horizontal). 
Two  times,  one  for  entry  into  the  beam,  (Tl) ,  and  one  for  exit 
from  the  beam,  (T2),  were  computed  for  each  target.  Thus,  the 
effective  antenna  pattern  was  "on”  over  a  fixed  angular  spread 
and  "off"  otherwise.  In  RSMOD  a  similar  computation  is  used  to 
ascertain  the  times  over  which  each  target  is  to  be  considered  in 
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a  broad  field  of  view.  Over  this  field  of  view,  the  antenna 
pattern  (presented  as  a  weighting  function) ,  target  cross-section 
(aspect  angle  and  depression  angle  dependent) ,  and  noise  level 
are  used  to  compute  an  amplitude  of  the  return  signal.  Transmitter 
and  receiver  gain  are  directly  related  to  the  antenna  pattern 
and  the  maximum  gain  of  the  antenna. 

Intermediate  to  the  receiver's  first  amplification  stage 
(receiver  gain)  and  the  video  recorder,  there  are  various  1  imitating 
factors:  noise  limiting  and  AGC  (automatic  gain  control)  are  two 
examples.  RF  limiting  will  act  with  no  time  delay.  When  limiting 
occurs,  the  limiting  level  will  be  a  function  of  the  total  power 
level  (signal  plus  background  noise  level),  limiting  levels  and 
nominal  amplifier  gain.  The  computed  gain,  under  limiting,  that 
is  applied  to  each  signal,  to  noise  level,  and  to  clutter  will 
preserve  S/N  ratios.  Similar  functions  or  possibly  tabulated 
weighting  functions  can  be  used  to  relate  signal  attenuation  through 
out  the  radar  processor  chain.  Computation  of  an  AGC  requires 
integration  of  returns  (at  the  video  stage)  from  several  targets. 

In  Fig.  8,  example  amplitudes  of  target  returns  from  a  number 
of  targets  are  plotted  relative  to  a  time  scale.  To  adjust  or 
set  the  AGC  from  this  data,  it  would  be  necessary  to  compute 
a  running  average  (integrate)  over  these  returns  under  the  in¬ 
fluence  of  a  time  constant  (AGC  lags  the  present  signal  strength) 
and  a  corresponding  delay  which  determines  the  weighting  function 
for  decay  of  olu  information.  Due  to  the  added  complexity  of  the 
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Figure  8.  Amplitude  of  Target  Return  vs.  Time 


AGC  correction,  we  have  not  included  this  consideration  in  RSMOD. 

Antenna  stabilization  errors  are  an  input  from  the  aircraft 
simulation.  Without  these  errors,  the  target  cross-section  profile, 
the  antenna  pattern,  the  instantaneous  squint  and  depression  angles, 
and  range  and  range  rate  to  the  target,  describe  an  expected  history 
of  amplitude,  video  frequency  an  doppler  frequency.  With  antenna 
pointing  errors: (a)  the  antenna  pattern  will  be  biased  from  and/ 
or  will  oscillate  about  the  expected  pattern  to  create  errors  in 
observed  amplitudes;  (b)  the  clutterlock  function  will  supply  either 
a  low  or  a  high  offset  frequency,  since  the  doppler  returns  from 
the  target  field  observed  will  be  shifted  away  from  the  desired 
doppler. 

MSTRIP  contains  a  time  slice  which  is  arbitrarily  based  on 
dividing  the  time  history  for  each  target  into  an  equal  number 
of  slices.  That  is,  if  a  target  is  in  the  beam  for  1.5  sec.  and 
the  number  of  slices  is  11,  then  a  target  history  will  be  generated 
in  each  0.15  sec.  beginning  with  T1  ■  tQ  and  going  to  T2  ■  tQ  + 

1.5.  In  the  updated  RSMOD*  PRF  timing  is  used.  Since  the  limiter 
functions  require  a  parallel  look  at  targets,  background  clutter 
and  noise,  common  time  points  are  required.  Time  slices  are 
taken  as  equally  spaced  (integral  multiples  of  1/PRF)  at  the 
beginning  and  then  initial  times  for  targets  are  adjusted  to 

*As  a  consequence  of  the  evolving  form  of  RSMOD,  it  is  of  some 
interest  to  note  that  MSTRIP  is  embedded  in  RSMOD.  Consequently, 
MSTRIP  and  other  simulation  forms,  intermediate  from  MSTRIP  to 
RSMOD,  can  be  operated  from  the  same  coded  program. 
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correspond  with  l’RR.  At  first  g Unco,  this  may  seem  to  bo  a  concern 
over  inconsequential*;  but,  it  is  a  reasonable  concern  since, 
doppier  spatial  frequency  In  the  range  direction  is  dependent 
oh  chirp  frequency  and  range  scale  factors.  The  doppier  is  also 
range  rate  dependent  and  will  depend  on  motion  compensation  which 
will  be  time  dependent  since  each  record  in  Che  range  direction 
begins  at  a  range  which  is  a  function  of  the  commanded  minimum 
range  and  is  also  a  function  of  the  recorder  sample,  rate. 

As  an  interaction  between  aircraft-target  dynamics  and  the 
radar  frequencies,  the  simulation  model  will  compute  time  his¬ 
tories  of  uncpmpehsa ted  doppier  for  each,  target  center,  Here¬ 
tofore,  these  dopplers  have  been  computed  as  a  function  of  the  radar, 
frequency;  a  single  fixed  frequency,  Mow,,  however,  with  the  addi¬ 
tion  of  chirp  frequencies,  t«a  doppier  is  dependent  on  the  chirp. 
R8M01)  computes  dopplers  as  a  function  of  chirp  frequencies. 

Doppler  at  the  leading  edge  of  the  chirp  pulse  is  based  on 


J(£c  + 


R,  and 


2(f  +  f  )  • 
o  n  R 


for  the  trailing  edge;  whore  t'c  »  center  frequency,  f^  *  chirp 
frequency  loading  edge,  C  »  chirp  frequency  trailing  edge, 

R  »  instantaneous  range  rate  to  target  in  question,  and  C  ■ 
velocity  of  propagation. 

6.2. 4  Motion  Compensation 

Inputs  to  the  computation  of  motion  compensation  have  been; 


nominal  altitude,  nominal  sweep-trigger  delay,  nominal  sweep  time, 
a  doppler-correction  function,  trigger  delay  function,  and  film- 
speed  function.  We  have  investigated  the  possibility  of  adding  a 
reference  path  as  an  explicit  input.  Through  the  use  of  AIRMOD 
this  capability  has  been  added  to  RSMOD.  Hence  aircraft  maneuver 
considerations  can  be  allowed.  In  the  MSTRIP  simulation,  a  single 
straight  line  was  assumed  as  the  reference  path.  If  AIRMOD  inputs 
are  used,  all  motion  compensation  can  be  brought  over  as  tabular 
entries.  Otherwise,  the  following  internal  computations  can  supply 
motion  compensation. 

A  phase-shifter  rate  (dopper  correction)  is  computed  as  a 
function  of  aircraft  motion  which  is  used  to  correct  the  com¬ 
puted  uncompensated  dopplers.  This  rate  is  obtained  by  choosing 
one  of  several  possible  motion-compensation  rules.  Currently, 
we  have  ten  individ  si  motion  compensation  rules  available  for 
examination.  Each  set  of  rules  specifies  the  functions  to  be  used 
to  compute  doppler  frequency  correction,  azimuth  sampling  rate, 
and  range-delay  correction. 

In  the  use  of  pure  broadside  (unsquinted)  strip-mapping, 
questions  of  proper  motion-compensation  rules  do  not  arise  (except 
for  the  relatively  straightforward  question  of  single-range  com¬ 
pensation  versus  multi-range  or  continuous-range  compensation, 
for  which  most  of  the  trade-offs  are  known).  For  squinted  (non¬ 
broadside)  strip-map  operation,  however,  determination  of  the 
best  signal  mix  for  motion  compensation  is  not  at  all  clear-cut, 
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except  for  the  case  of  a  single  point  target  at  zero  depression 
angle.  Moreover,  even  the  one-degree  squint  angle  resulting  from 
the  use  of  azimuth  (doppler)  offset  in  "broadside"  mapping  can  be 
appreciably  non-zero  at  steep  depression  angles.  Consider  the 
forward-squint  case,  with  single-range  motion  compensation. 

The  sketch,  Fig.  9,  shows  velocity  components  in  four 
directions:  beam  velocity  (V0)  down  the  squinted  radar  beam 

D 

to  swath  center;  cross-beam  velocity  (V,,,.)  normal  to  the  radar  beam; 
along-track  velocity  (VT)  in  the  direction  of  the  reference  path; 
and  cross-track  velocity  (V^)  normal  to  reference  path.  In 
the  broadside  case,  V^,  and  are  identical,  and  this  quantity 
is  used  to  control  azimuth  sampling  (thus  locking  actual  dis¬ 
placement  to  ground  displacement,  to  eliminate  along-track 
geometric  distortion  in  the  map);  similarly,  V  and  V  _  are 
identical,  and  this  is  the  quantity  required  for  frequency- 
correction  of  the  doppler  returns  (to  eliminate  cross-track 
distortion  and  smear  in  processing).  But  in  the  squinted  case, 
which  quantity,  Vg  or  V^,  should  be  used  for  doppler  correction? 
And  if  Vn  is  chosen,  should  not  the  orthogonal  velocity,  V^,,, 
be  used  instead  of  to  control  azimuth  sampling?  (If  so, 
azimuth  sampling  is  no  longer  slaved  to  ground-speed).  Which¬ 
ever  combination  is  chosen,  what  happens  to  the  image  quality 
cf  targets  appreciably  distant  from  the  swath  center? 

The  simplified  problem  of  a  single  point  target  at  swath 
center,  and  zero  depression  angle,  has  been  solved  analytically — 


the  answer  in  this  case  is  that,  for  best  azimuth  resolution, 
should  be  used  for  doppler  correction  and  V__  for  control  of 
azimuth  sampling.  But  the  more  realistic  cases  of  non-zero  swath 
width  and/or  non-zero  depression  angle  are  not  easily  analyzable — 
it  was  specifically  for  the  purpose  of  studying  these  cases  that 
the  SAR  simulation  model  was  originally  developed.  A  variety 
of  motion-compensation  rules  can  be  entered  as  alternate  subroutines. 
(For  examp] e,  MOCOMP  1  might  use  Vg  for  doppler  correction,  V^g 
to  drive  azimuth  sampling,  and  cross-track  displacement  for 
range-delay  correction;  MOCOMP  2  might  use  V  for  doppler,  V_ 

D  JL 

for  azimuth  sampling,  and  beam-direction  displacement  for  range 
delay.)  The  same  target  array  can  be  mapped  (with  the  same 
aircraft  dynamics)  alternately  using  each  of  the  MOCOMP  methods, 
and  the  resulting  processed  imagery  compared  for  resolution  and 
geometric  fidelity. 

The  motion-compensation  technique  study,  described  above, 
is  not  yet  complete.  Preliminary  results,  however,  seem  to 
indirate  that,  for  range  swaths  exceeding  10%  of  swath-center  range, 
the  best  combination  of  signals  for  single-range  compensation  is 
that  using  Vg  for  doppler  correction,  for  azimuth  sampling, 
and  displacement  in  the  beam  direction  for  range-delay  correction. 
This  compensation  method  degrades  resolution  somewhat  at  swath 
center,  but  resolution  at  swath  edges  is  better  preserved. 

Heretofore,  these  functions  have  been  computed  x’elative  to 
perceived  or  measured  motion  of  the  aircraft.  Mc'zion  measurements 
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and  generation  of  motion  compensation  did  not  have  error  sources 
as  inputs.  The  only  errors  that  existed  were  involved  with  the 
motion  compensation  rules  and  the  amount  of  uncompensated  motion 
that  ended  up  on  the  stored  returns,  plus  an  error  which  was  allowed 
as  the  difference  between  actual  forward  motion  and  measured  forward 
motion  of  the  aircraft.  In  the  updated  model  these  computations 
will  allow  input  of  RMS  errors  (deterministic,  not  random).  Thus, 
motion  may  be  measured  with  errors  and  motion  compensation  can  be 
computed  with  error.  For  these  reasons,  as  with  the  actual  systems, 
the  accumulation  of  motion  compensation  errors  and  residual  uncom¬ 
pensated  motion  can  cause  final  images  to  be  degraded  below  a 
minimum  acceptable  level  which  is  to  be  determined. 

The  function  of  the  clutterlock  is  to  generate  a  correction 
(error  signal)  to  bring  the  doppler  observed  after  motion  compen¬ 
sation,  as  measured  against  the  target  background,  into  coincidence 
with  a  desired  band  of  doppler  frequencies.  To  provide  this  result, 
a  clutterlock  can  be  designed  to  shift  all  doppler  frequencies, 
after  phase  shifter  compensation,  so  that  any  target  in  the  center 
range  position  at  the  aft  edge  of  the  antenna  beamwidth,  for 
example,  will  record  at  zero  frequency  or  will  record  at  a  prescribed 
minimum  frequency.  Other  choices  of  target  background  are  to 
average  over  a  band  of  possible  target  positions,  or  to  average 
over  several  bands  of  possible  target  ranges  and  azimuth  positions. 

A  residual  error  will  also  be  allowed  to  correspond  to  errors 
in  the  clutterlock  or  to  delays  in  making  the  corrections,  since 
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the  clutterlock  has  a  time  constant,  settling  time,  which  nominally 
is  fairly  long  (several  beamwidths).  Alternatively,  clutterlock 
signals  (with  errors)  can  be  inserted  from  the  AIRMOD  file. 

Since  we  are  providing  a  representation  of  chirp  and  are 
allowing  doppler  returns  to  be  chirp-dependent,  we  have  included 
the  effect  of  a  range  dependent  phase  shifter  on  chirp  frequencies. 
Hence,  chirp  warp  is  represented  by  a  function  generator  in  the 
simulation  model. 

6.2.5  Signal  Generation  and  Data  Storage 

Generation  of  signal  histories  depends  on  target,  aircraft, 
radar,  and  motion  compensation  inputs.  Also,  we  observe  that 
a  number  of  logic  variables  are  brought  into  play.  Some 
examples  are: 

a.  If  LTABMD  *  .TRUE,  then  changes  in  doppler  and  chirp 
frequencies  due  to  motion  compensation  may  be  obtained  from 
AIRMOD  tables  depending  on  whether  or  not  LPSRT,  LFLMSP,  LCUT 
and  LTDEL  -  representing  phase  shifter  input,  film  speed  input, 
clutterlock  input,  and  recorder  trigger  delay,  respectively  - 
are  TRUE  or  FALSE.  True  means  'use  AIRMOD'  and  false  means 
'use  internally  generated  function.' 

b.  If  LANG  *  TRUE,  target  motion  inputs  are  tabular  (e.g., 
data  logs  from  field  experiment).  Otherwise  function  input 

is  used. 

c.  If  LCHIRP  =  FALSE,  one  dimensional  signal  history  is 
generated  and  stored  corresponding  to  doppler  returns  only. 
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Otherwise,  two  dimensional  signal  history  chirp  (video)  in  range 
and  doppler  in  azimuth  -  is  generated  and  stored.  If  LMAINI  is  also 
turned  on,  equals  TRUE,  then  appropriate  binary  records  are 
written  into  a  file  for  future  use  in  computing  impulse  responses. 

A  total  of  fifteen  logical  (control)  parameters  are  initial¬ 
ized  by  read  in  of  the  first  data  card.  For  the  purpose  of  print¬ 
ing  out  intermediate  results  during  debugging  operations,  there 
are  an  additional  nineteen  debug  controls  read-in  on  the  second 
card.  Control  cards  are  most  important  with  respect  to  actual 
use  of  the  simulation  program.  Further  discussion  of  the  con¬ 
trol  or  logical  inputs  is  subjugated,  for  the  most  part,  to  the 
details  of  program  running  that  are  presented  in  the  Appendix 
on  RSMOD  programs. 

A  recording  process  to  produce  a  digitized  signal  is  em¬ 
ployed  in  MSTRIP,  which  computes  azimuth  spatial  frequency  from 
compensated  doppler  and  commanded  azimuth  sampling  rate.  These 
doppler  histories  are  attached  to  storage  locations  in  azimuth 
and  in  range  by  a  computed  data  position.  Data  positic*  is  a 
function  of  a  commanded  range-delay  trigger,  range  sampling  rate, 
and  azimuth  sampling  rate.  We  have  added  to  these  histories  a 
corresponding  spatial  frequency  in  the  range  direction  that 
results  from  the  received  chirp  signals.  Since  the  chirp  fre¬ 
quencies,  added  to  the  carrier  frequency,  produce  a  variation 
in  doppler  across  the  chirp,  the  new  recorder  attaches  dopplers 


for  each  chirp  pulse  that  are  chirp-frequency  dependent.  Thus 
the  effect  of  chirp  on  azimuth  (doppler)  frequency  is  represented. 
We  have  a  simulated  data-storage,  then,  which  records  video  fre¬ 
quencies  in  the  range  direction  and  doppler  frequencies  in  the 
azimuth  direction.  Spatial  frequencies  in  the  range  (video) 
direction  varies  with  range  sampling  rate.  Sampling  rates  are 
allowed  to  have  slowly-varying  and  impulse-type  nonlinearities. 

Another  important  feature  of  RSMOD  is  target  amplitude 
(signal  strength  of  target  returns).  Gain  and  limiting  functions 
were  mentioned  earlier  under  the  radar  sub-section.  After  the 
signal  strength  at  the  receiver  front-end  has  been  attenuated 
by  noise  limiting,  and  other  system  attenuations,  we  take  the 
resulting  signals  and  record  a  corresponding  intensity  for  each 
data  point.  There  is,  therefore,  a  history  of  recorded  signal 
and  noise  amplitudes  as  functions  of  video  amplifier  bandwidth 
(and  CRT  spot  size  for  optical  recording)  in  range  and  in  azimuth 
(recorder  MTF  characteristics). 

Finally,  the  recorded  signal  histories  can  provide  a  two- 
dimensional  "signal-film"  output  in  a  format  for  impulse-response 
calculations. 

6.2.6  Processor 

A  processor  exists  in  RSMOD  (from  MSTRIP)  which  accepts 
the  following  input  parameters,  which  are  characteristics  of 
the  processing: 
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a.  Focus-target  selection 

b.  Best-focus  criterion  (minimum  azimuth  smear,  minimum 
range  smear,  minimum  real  mean  smear) 

c.  Squint-angle-correction  criterion 

d.  Depression-correction  criterion 

e.  Fixed-frequency  gate  option. 

These  parameters  enter  into  the  azimuth  offset  processing  of 
doppler  histories  (azimuth  focus).  If  chirp  (video)  signals  are 
to  be  focused  a  new  processing  component  called  RNGPRC  for  range 
processor  is  invoked.  In  which  case,  individual  video  signals 
for  a  target  are  focused  (compressed)  to  form  a  one-dimensional 
signal  history.  This  range  compressed  history  is  identical  in 
form  to  the  azimuth  histories  used  in  MSTRIP.  Therefore,  the 
second  stage  of  processing  for  RSMOD,  is  identical  with  MSTRIP 
processing.  Before  processing  any  of  the  doppler  histories, 
the  processor  is  primed  with  a  fixed  band  of  spatial-frequency- 
response  limitations.  Then  four  optimization  procedures  are 
available.  Three  of  these  are  computed  automatically  but  are 
dependent  on  the  selection  of  a  focus  target  (one  of  the  input 
targets)  and  the  selection  of  a  best-focus  criterion.  In  these 
cases,  the  focus  target  and  its  location  in  range  and  azimuth 
(its  doppler  and,  eventually,  video  frequency  history)  becomes 
the  controlling  factor  for  computing  smear  and  location  of  other 
targets.  Targets  are  excluded  or  included  according  to  the 
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overlap  of  signal  histories  (in  the  frequency  domain)  with  the 
focus  target.  A  target  can  be  focused  if  and  only  if  a  minimum 
history  criterion  is  met.  Presently  this  criterion  says  that  at 
least  14%  (3  points  out  of  21)  of  a  target's  doppler  history  must 
overlap  with  the  optimized  frequency  band  and/or  the  fixed  fre¬ 
quency  band  defined  for  first  looks  at  the  focus  target.  If  the 
focus  target's  doppler  bandwidth  does  not  meet  this  criterion, 
no  targets  are  focused  and  a  comment  is  printed  to  the  effect 
that  the  focus  target  has  no  useable  doppler  history.  The  focus 
target  signal  history  is  printed  out,  then,  and  the  program 
returns  to  the  next  computer  run. 

A  fourth  method  of  optimization  is  allowed — namely  a  fixed 
spatial  frequency  band  can  be  stipulated  as  before,  but  this 
frequency  band  can  be  retained  for  checking  and  focusing  the 
signals  for  each  target. 

If  only  MSTRIP  is  used,  the  processor  section  focuses  in 
azimutn  only.  It  computes  azimuth  and  range  smear  based  on  using 
data  already  compressed  (focused)  in  the  range  direction.  Since 
chirp  is  already  compressed,  (short  pulse  system)  nonlinearity  of 
chirp  implications  on  range  resolution  is  fed  in  as  a  fixed  smear, 
RRMIN,  due  to  chirp.  This  RRMIN  input  is  then  averaged  (RSS) 
with  nonlinearities  due  to  "range  walk"  to  obtain  a  range  reso¬ 
lution.  For  the  expanded  version  of  RSMOD  ,  compression  of  chirp 
and  the  accumulated  effects  of  errors  in  generating,  transmitting, 
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receiving,  and  recording  chirp  signals  are  computed  during  each 
simulation  run.  In  which  case,  RRMIN  becomes  a  computed  number  for 
each  range  compression.  Since  range  compression  can  now  produce 
the  doppler  history  for  the  second  stage,  there  is  a  corresponding 
contribution  to  doppler  noise  or  azimuth  smear. 

Outputs  from  the  processor  are: 

a.  A  list  of  all  input  parameters 

b.  Range  and  azimuth  image  position  for  each  target 

c.  Range,  azimuth,  and  average  (RMS)  smear  for  each 
target. 

d.  Geometric  distortion  between  targets,  with  and  without 
corrections. 

e.  Printplot  of  output  image  with  scale  factors. 

6.2.7  Displays 

Output  of  video  histories,  doppler  histories,  processor 
results,  plotted  image  positions,  input  parameter  lists,  binary 
inputs  for  impulse  response,  and  various  intermediate  parameter 
values  are  controlled  by  integer  and  logic  switches.  The  switches 
are  input  parameters.  Logic  structure  in  the  four  sections  of 
the  RSMOD  program  also  control  display  of  information  in  terms 
of  existence  of  a  particular  kind  of  record.  For  example,  if 
LCHIRP  is  off,  no  video  history  is  computed.  As  a  consequence 
of  this  switch  setting  the  program  never  enters  the  write  out 
section  in  which  video  histories  are  displayed.  Hence,  the 


writeout  switch  for  2D  video-doppler  display  could  be  on  but 
the  output  of  video  would  not  be  attempted  unless  LCHIRP  is 
also  on. 

6.3  Signal  Generation 

Signal  generation  as  a  module  in  the  simulation  program  consists 
of  a  main  module  (SIGGEN)  which  calls  sub-modules  and  stores  signal 
history  according  to  logic  controls.  Choices  of  limiters,  time 
base,  limiter  computations,  and  one  or  two  dimensional  histories 
are  represented  in  these  controls.  During  the  signal  generation 
stage,  four  important  sub-modules  are  involved.  These  are: 

a.  TIMER  -  a  calculation  of  target-radar  Interaction  times 

b.  POINT  -  a  collection  of  algorithms  which  supplies  trans¬ 
lations,  rotations,  and  vector  calculations  of  target 
points  relative  to  the  phase  center  of  the  radar 

c.  MOCOMP  -  a  group  of  computations  which  provide  motion 
compensation  values  for  the  determination  of  signal 
history 

d.  TABLE  -  a  module  which  generates  and  tabulates  signal 
history  by  using  results  from  TIMER  and  calculations 
from  POINT  and  MOCOMP. 

In  TABLE,  signal  amplitudes  are  calculated  to  the  receiver  front- 
end  and  chirp  images  are  developed  for  video  returns.  TABLE 
works  on  an  individual  target  basis  while  SIGGEN  associates  these 
individual  scalar  and  vector  results  with  time,  and  stores  values 


-46- 


in  appropriate  arrays  indexed  by  target,  by  time  slice,  by  range 
slot,  and  by  chirp  sample. 

Figure  10  illustrates  the  order  of  computation  and  the 
options  available.  These  options  are  determined  by  logical  switches 
which  are:  LAMP  (for  amplitude  computation),  LNEWT  (for  timing 
sequence),  LCHIRP  (either  azimuth  only  signals  or  Range  and  Azimuth 
Signals),  LLIMIT  (determines  limiting  for  RF) ,  and  LVDAMP  (corre¬ 
sponding  to  video  amplification).  Other  options  appear  in  TABLE, 
M0C0MP,  and  POINT  relative  to  tabular  inputs  vs.  function  input 
and  choice  of  procedure.  We  will  now  discuss  the  options  and  re¬ 
sulting  computations  in  the  context  of  the  various  subroutines 
and  function  routines  that  are  used. 

6.3.1  TIMER 

TIMER  is  the  subroutine  in  which  the  entry  and  exit  time 
for  each  target  is  computed.  Either  a  scanning  beam  or  a  scanning 
beam  and  a  looking  beam  may  be  specified.  The  algorithm  for  the 
determination  of  these  times  is  based  on  squint  angles,  0^>  G^, 

@2,  and  38  shown  in  Fig.  11.  For  purposes  of  illustration, 
we  let  the  target  point  T  walk  through  the  beam  along  the  dotted 
line  as  shown  and  assume  the  velocity  of  T  to  be  the  real  velocity 
of  the  target  minus  the  aircraft  ground  speed.  @T  the  initial 
squint  angle  to  T  is  computed  by  using  POINT.  Then,  if  the  time 
to  0^,  for  example,  is  desired,  an  estimated  time  to  close  ' 
computed  from  ground  range  to  target,  aircraft  ground  speed,  and 


-47- 


♦ _ 

CALL  TABLE  TO 
PRODUCE  A  POINT 
AZIMUTH  OF  HISTORY 
FOR  THIS  TARGET 


Y 


IF  LIMITING  FUNCTION  THEN 
ACCUMULATE  TARGET  AMPLITUDE 
AND  ACCUMULATE  BACKGROUND 
ON  REQUEST 


RECORD  ID  (AZIMUTH) 
HISTORY 


WRITE  OUT  COMPLETE 
20  SIGNAL  HISTORY 
IF  DESIRED 


! 

i 


) 


Figure  10  (b)  SIGGEN  Generalized  Flow  Diagram  (continued). 
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the  angular  separation.  This  process  Is  iterated  three  or  four 
times  until  the  computed  9  at  some  location  like  X  is  near  enough 
to  0^.  We  observe  that  aircraft  and  target  velocities  vary  with 
time  and  hence  an  iterative  procedure  is  required. 

6.3.2  TABLE 

The  subroutine  TABLE  provides  control  and  computations  to 
generate  signal  histories  on  a  per  target  and  a  per  time  slice 
basis.  Signal  history  in  chirp  and  doppler  (range  and  azimuth) 
is  computed  from  transmitted  chirp  frequencies,  beam  velocity, 
and  applicable  motion  compensation.  It  is  these  chirp  and  doppler 
histories,  considered  it;  a  spatial  orientation  that  are  compressed 
(focused)  to  obtain  final  images  and  resolutions.  In  addition  to 
these  spatial  frequencies,  amplitude  of  target  returns,  background 
clutter,  and  noise  figures  are  also  computed.  Since  the  subroutine 
TABLE  involves  so  many  interactive  elements  it  is  desirable  to  show 
a  generalized  flow  diagram  at  this  hierarchical  level  also.  In 
Fig.  12,  it  is  observed  that  subroutines  POINT,  MOCOMP,  and 
CHIRPX  are  used  along  with  a  function  named  AMPR. 

6.3.3  POINT 

Subroutine  POINT  contains  aircraft  motion  and  target  motion 
functions.  These  are  quadratic  functions  of  time  with  sinusoid 
and  step  function  terms.  Velocity  equations  and  angle  generating 
function  for  target  roll,  pitch  and  yaw  are  included.  Alternatively 
within  this  subroutine,  calls  on  AMOD  and  SHIPMO  subroutines 
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Figure  12.  (a)  Subroutine  Table 


Figure  12  (b)  Subroutine  Table  (continued). 


-53- 


provide  more  realistic  aircraft  and  target  motion.  As  part  of 
this  alternative,  one  may  specify  table  look-up,  linear  inter¬ 
polation,  or  quadratic  interpolation.  It  is  in  POINT  and  in 
MOCOMP  where  the  aircraft  simulator  (or  aircraft  logs)  are 
used  extensively. 

6.3.4  MOCOMP 

Subroutine  MOCOMP  contains  the  functions  and  table  look¬ 
ups  which  perform  motion  compensation.  Motion  compensation 
was  discussed  earlier  in  terms  of  input  parameters  and  relative 
to  parametric  studies  aimed  at  finding  best  motion  compensation. 
With  the  addition  of  AIRMOD  inputs  to  RSMOD,  we  can  now  provide 
a  better  input  of  error  sources  as  well  as  different  forms  of 
motion  compensation.  Heretofore,  (MSTRIP) ,  we  studied  motion 
compensation  relative  to  aircraft  motion  or  velocity  errors 
versus  different  rules  for  compensation.  Now  error  sources 
(RMS)  can  be  introduced  realistically  into  AIRMOD  and  relayed 
via  phase  shifter,  range  delay  time,  azimuth  sample  rate,  clutter- 
lock  and  antenna  pointing  errors.  Moreover,  the  effect  of  INS 
(Inertial  Navigation  System)  errors  is  introduced  via  a  real  INS 
representation  in  AIRMOD.  Since  the  error  generators  are  fully 
controlled  by  the  AIRMOD  simulator,  all  errors  with  respect  to 
time,  as  well  as  the  effects  of  errors,  can  be  an  output  for 
subsequent  analyses. 

6.1.5  AMPR 

Computing  amplitude  of  target  returns  is  accomplished  through 

-54- 


v> 'X ---y-rygT y' 


--  -= a 


Che  use  of  a  number  of  subroutines  and  function  subroutines. 

To  compute  amplitude,  in  TABLE,  it  is  only  necessary  to  assign 
AMPR,  a  function  of  target  number,  time,  unit  range  vector, 
slant  range  to  target,  and  an  attenuation  factor,  to  an  appro¬ 
priate  variable.  In  the  coded  program,  the  Fortran  Statement 
is 

IF  (LAMP)  AMP  =  AMPR(NTP,  TT,  UR,  RA,  ATEN) 

Hidden  from  view,  in  AMPR,  are:  the  antenna  and  target  pattern 
functions;  antenna  and  target  pointing  functions;  the  radar  power 
equation;  tables  of  antenna  patterns;  tables  of  target  cross- 
section,  pointing  angles,  and  lobe  patterns;  a  sine  function 
generator.  Detailed  equations  for  the  AMPR  function  are  con¬ 
tained  in  the  coded  programs  that  are  listed  in  the  RSMOD  appendix. 
Here  again  AMOD  tables  from  AIRMOD  provide  an  important  input 
by  relating  antenna  pointing  angles  relative  to  aircraft  frame 
with  respect  to  time.  These  angles  combined  with  aircraft  roll, 
pitch  and  yaw  (from  AIRMOD)  determine  the  instantaneous  squint, 
depression  and  pattern  orientation  for  the  antenna.  There  are  times 
when  it  may  be  desirable  to  back-off  from  AIRMOD  realism.  One 
may  wish  to  study  some  special  effect  or,  during  a  phase  of  program 
development  and/or  debugging,  static  antenna  patterns  may  be 
desirable.  For  this  purpose,  idealized  routines  are  provided, 
via  switch  controls,  to  give  aircraft  and  antenna  rotations, 
either  fixed  or  varied  with  time. 
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Functions  and  routines  which  support  AMPR  are  designed  to 
be  easily  reprogrammable.  For  example,  the  antenna  pattern,  pre¬ 
sently  used,  is  the  product  of  a  vertical  and  a  horizontal  pattern. 
To  compute  antenna  gain  along  a  particular  direction,  these  patterns 
are  sampled  and  the  results  multiplied  together  (added  in  DB) . 

A  vertical  antenna  pattern  is  shown  in  Fig.  13.  Depicted  is  a 

2 

portion  of  a  calibration  curve  obtained  from  a  esc  antenna 
pattern.  For  simulation  purposes  the  peak  power  point  is  shown 
at  0°.  Positive  depression  (down)  is  plotted  to  the  left  of 
zero.  In  use  the  pointing  angle  of  the  antenna  a  pseudo  bore 
sight,  corresponds  to  the  zero-angle  gain  (max.  gain)  shown. 

A  maximum  gain  of  30  db  above  isotropic  is  representative. 


2  /sin  x2\ 

Using  a  sine  x  =  I  function  for  the  horizontal  pattern 

produces  a  gain  pattern  similar  to  Fig  14.  If  0  is  the  hori- 

7TD 

zontal  aspect  angle,  then  x  =  —  sin  0,  where  A  is  wavelength 
and  D  is  antenna  aperture.  In  Fig.  14.  If  0  is  the  hori- 
.espect  to  sin  0.  The  half-power  beam  width  is  singled  out  as 
.88  For  convenience  in  operating  the  model  and  since  —  ■'s 
often  used  as  an  approximation  for  the  rectangular  equivalent 
half-power  beamwidth,  the  simulation  mechanizes  sine  functions 
with  the  half-power  beams  as  an  input.  As  a  consequence  of 
choosing  this  input,  antenna  patterns  and  cross-section  of  tar¬ 
gets  can  be  idealized  by  the  same  kind  of  sine  functions.  0  = 

1°  would  be  a  nominal  beamwidth  for  the  radar  antenna  at  X-band 
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and  0  *  10°  (or  30°)  might  be  representative  for  a  target  cross- 
section  pattern.  (It  is  not  desirable,  at  this  time,  to  entertain 
any  kind  of  rigorous  representation  of  target  signatures  in  a 
systems  model.  Too  much  rigor  in  an  otherwise,  relatively  simple 
representation  would  not  be  beneficial). 

Before  leaving  our  discussion  of  antenna  patterns,  one 
additional  feature  can  be  used  to  vary  the  horizontal  structure. 
Coded  into  the  present  formulation  we  admit  a  weighting  function 
in  tabular  form  which  can  be  used  to  describe  antenna  patterns 
that  are  asymetrical  in  horizontal  structure. 

Inputs,  for  each  target,  describe  beamwidths  for  horizontal 
and  vertical  lobing  of  cross-section,  maximum  cross-section,  and 
the  roll-pitch-yaw  orientations  of  the  cross-section  pattern. 

Each  point  is  coordinatized  with  respect  to  the  target  frame, 
which  may  be  angularly  oriented  with  time.  Henca,  aspect  angle 
dependence  as  seen  by  the  radar,  may  vary  with  time. 

The  following,  independely  coded,  functions  with  outputs 
noted  are  used  to  compute  AMPR: 

a.  Either  AIRANG  (AIRMOD)  or  ANGAIR  (function)  -  angular 
orientation  of  the  aircraft 

b.  Either  ANTPNT  (AIRMOD)  or  ANGANT  (function)  -  squint 
and  depression  angle  of  the  antenna  relative  to  the 
air  frame 

c.  ANPNT  -  rotation  of  antenna  axes  to  air  frame 
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d.  ANTDEP  -  verticle  antenna  pattern 

e.  ANTSQT  -  horizontal  antenna  pattern 

f.  TGTSIG  -  target  cross-section 

g.  RADRDB  -  Anplitude  of  a  target  return  in  DB. 

6.3.6  CHIRPX 

In  CHIRPX,  the  transmitted  chirp  pulse  is  generated.  A 
beginning  chirp  fxequency,  a  terminal  chirp  frequency,  a  pulse 
width,  and  a  sample  rate  are  inputs.  From  these  inputs,  NCH1RP 
(the  number  of  chirp  samples  for  simulation)  linear  samples  of 

the  chirped  pulse  are  formed.  From  either  a  table  (read  in)  or 

# 

a  function  generator  (coefficients  read  in);  NCHIRP  samples  of 
a  non-linear  function  are  added  to  each  linear  sample  t)oint  to 
form  the  non-linear  pulse  that  is  transmitted.  At  this  time,  the 
only  non-linearities  which  can  appear  in  the  chirp  as  recorded 
ares  transmitter  generated  and  those  which  result  from  the  geo¬ 
metry  and  motion  compensation.  That  is,  we  have  not  added  phase 
shifts  due  to  the  radar  receiver  or  any  other  function  downstream 
to  the  recorder.  Obviously  some  phase  shifting  occurs  but  it  is 
assumed  to  be  of  second  order  importance  to  those  represented  in 
the  model.  We  have  already  discussed,  phasing  shifting  and  chirp 
warp  in  the  input  section  above. 

6.3.7  LIMIT 

After  table  generation  is  accomplished  over  all  targets  for 
a  given  time  slice,  the  gain  and  limiter  functions  are  performed. 
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Either  RF  amplifier  or  video  amplifier  gains  can  be  combined 
with  limiter  functions  by  the  following  representation.  Consider 
an  amplifier  gain  function  as  shown  in  Fig.  15.  We  have  shown 
a  piece-wise  linear  function  with  two  pieces.  The  first  linear 
piece,  corresponding  to  Gain  1,  relates  input  power  to  output 
power  for  non-limiting  input  power.  If  the  total  signal  from 
noise,  targets  and  background  exceed  X^m»  the  power  input  level 
over  a  pulse  width  at  which  limiting  begins,  then  the  following 
considerations  hold.  First,  after  limiting  begins,  signal-to- 
noise  (S/N)  input  ratios  must  be  equal  to  S/N  output.  That  is, 
a  good  approximation  to  limiting  in  an  RF-amplifier  or  video¬ 
amplifier  is  obtained  if  a  linear  model  is  used  in  whici:  S/N 
is  preserved.  Secondly,  the  gain  after  limiting  for  the  total 
signal  is  to  be  GAIN2,  while  GAIN1  holds  prior  to  limiting.  These 
two  conditions  are  sufficient  to  establish  that: 


y  -  GAIN1  *  x,  X„  „  .  <  X..  , 
total  —  lim 

is  the  equation  of  the  first  linear  piece 


(ii) 

(iii) 


y..  -  GAIN1  *  X., 

7 lim  lim 

Y„  .  .  -  GAIN2  *X.+Y,X>X,. 
total  total  o  lim, 


with 
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which  implies  that  for  all  X,  when  X  >  X.^^, 

y  •  GA1N3  *  X  -  1GA1M2  *  Xtotal  +  Xli„  (GAIN1  ~  GA1N2)1  *X 

^total 

and 

GAIN3  «  GATN2  *  Xtotal  +  Xlim  (GAIN1  ~  GAIN2> 

X 

total 

also 

Noise-out  -  (NOISE-IN)  *  GAIN3 
6.3.8  MTF's 

As  indicated  earlier,  MTF's  (Modulation  Transfer  Function), 
particularly  in  the  recording  process  for  an  optical  processor, 
are  important  contributors  to  recorded  amplitudes.  It  is  impor¬ 
tant  in  any  validation  of  this  simulation  model  that  emperical 
results  can  be  interleaved  with  simulation  results.  We  have 
chosen,  therefore,  to  model  an  optical  recorder.  In  only  these 
MTF's  and  in  the  linear  processor  of  the  next  section  do  we  find 
specialization  to  optical  processing.  All  other  representations 
including  the  impulse  response  calculations  are  not  restricted  in 
this  way.  By  disconnecting  these  MTF's  and  substituting  a  par¬ 
ticular  digital  processor  representation,  the  present  RSMOD 
should  modify  quite  readily  into  a  simulation  of  a  digital  processor. 

There  are  four  MTF's  which  have  been  singled  out  so  far. 

These  are  CRT-recorder ,  azimuth  and  range,  response  to  temporal 
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and  spatial  frequency,  lens  MTF  and  film  MTF.  Figures  16,  17, 
and  18  are  samples  of  calibration  curves  for  three  of  these 
MTF’s.  Since  response  curves  for  these  transfer  functions  are 
available,  we  can  use  tabulated  values  to  represent  these  MTF’s. 

The  basic  relations  for  spot  size  and  spatially  dependent 
transfer  functions  are  illustrated  in  Figure  19.  At  an  input 
spatial  frequency  (a  function  of  the  temporal  frequency  and 
sample  rate),  the  attenuation  begins  to  increase  as  a  quadratic 
function  of  temporal  frequency  and  continues  to  an  input  of 
(X^  <  X^) .  At  which  point  spot  size  and  spatial  frequency  becomes 
the  dominate  factor  in  determining  attenuation.  From  X£  to 
larger  X's,  the  attenuation  (db)  increases  almost  linearly  and  is 
representable  as  a  linear  function.  A  point  on  the  curve  at  X^ 
is  required  to  fit  the  quadratic  function.  Without  including  a 
specific  model,  dependent  on  spot  size,  we  can,  therefore,  include 
functions  which  are  described  by  inputs  X^,  ^3*  anc*  ^4  along 

with  the  level  of  attenuation  at  X^,  and  X^.  In  this  manner 
spot  size  is  included  implicitly.  A  logic  variable  in  each  of 
the  coded  transfer  functions  enables  switching  back  and  forth 
between  a  table  and  a  function. 

6.4  PROCESSOR 

6.4.1  Focusing  Signal  Histories 

Our  processor  for  signal  histories  in  optical  form  is  con¬ 
structed  through  the  use  of  a  linear  fit  procedure.  Consider 
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Figure  16.  Range  Response  MTF. 
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SPATIAL  FREQUENCY  (on  the  film)  IN  OPTICAL  CYCLES/mm 


Figure  19.  Recorder  Limiting  Function. 


che  print  out  of  signal  history  as  shown  in  Table  1  and  columns 
FSPAT  (FS)  and  z  which  denote  the  spatial  frequency  and  the  posi¬ 
tion  along  the  doppler  history  for  target  point  2.  For  FS  vs.  Z 
to  correspond  to  a  perfect  section  of  signal  film  (will  focus 
to  a  point  in  the  image  plane),  FS  must  be  a  linear  function  of  Z. 
The  sample  plot  with  circles  in  Fig.  20  corresponding  to  Table 
I,  is  nearly  linear.  This  corresponds  to  having  quadratic 
phase  in  a  recorded  signal  of  phase  history.  The  amount  of 
non-linearity,  deviation,  (DEVA) ,  of  the  actual  points  of  FS  to 
the  best  linear  fit  of  FS  vs.  Z  is  a  measure  of  how  much  azimuth 
smear  can  be  attributed  to  this  non-linearity.  Azimuth  smear 
SMA  from  DEVA  is  then  DEVA/ (slope  of  linear  fit)  which  provides 
the  contribution  of  FS  deviation  in  the  determination  of  G  which  is 
essentially  in  the  direction  of  Z.  The  second  plot  in  Fig.  6.16, 
shown  with  A's,  for  similar  geometry  but  with  a  lot  of  uncompen¬ 
sated  motion  is  more  illustrative  of  smear  as  predicted  by  a 
linear  fit.  Another  important  contributor  to  resolution  is  the 
amount  of  useable  signal  history.  Useable  is  defined  in  the  next 
sub-section.  If  the  entire  spatial  history  for  a  target  is 
useable,  then  the  effective  synthetic  aperture  (DEFF)  is  identical 
to  the  actual  aperture  (DACT) ,  where  DACT  ■  1/ (FSCALE*BEAM)  with 
FSCALE  *  2/A  and  BEAM  ■  A  is  wavelength  and  j  is  the  physical 
aperture  of  the  antenna.  D  might  be  5  to  10  feet.  Hence  for 
X-band  (A  *  .1  ft)  DACT  would  be  2.5  ft  to  5  ft.  As  is  well 
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Figure  20,  Sample  Plots  FS  vs 


knowu ,  0/2  is  the  theoretical  best  azimuth  resolution.  To  illus¬ 
trate  the  effect  of  useable  information  content,  let  21  points 
of  signal  history  represent  the  whole  synthetic  aperture.  If, 

the  largest  useable  history  were  10  points,  for  example,  then 
20 

is  the  degradation  factor  due  to  loss  of  information.  In 
other  words  the  effective  aperture  would  be 

DEFF  -  DACT*20/9 

or  DEFF  would  take  on  a  value  of  11  to  22  feet.  In  our  simula¬ 
tion,  this  contribution  to  smear  would  be  summed  (root  sum  of 
squares)  with  the  contribution  to  smear  given  by  the  SMA.  Thus 
azimuth  resolution  is  computed  as 

V(DEFF)2  +  (SMA)2 

A  similar  linear  fit  for  W  vs  G,  W  is  range  position  and  G 
is  azimuth  positicr ,  provides  a  smear  or  contribution  (DEVR)  to 
range  resolution  that  results  from  range  walk  and  uncompensated 
motion  of  aircraft  or  target.  The  other  major  contributor  to 
range  smear  is  chirp  non-linearity  (similar  to  FS  non-linearity). 
After  recording  chirp,  temporal  frequencies  relative  to  a  range 
sampler,  chirp  frequencies  appear  as  spatial  distributions  FSV. 
These  FSV  distributions  are  distributed  in  range.  Hence,  a  linear 
fit  of  FSV  vs.  WC,  where  WC  is  the  spatial  distribution  of  a 
chirp  pulse  in  range,  will  provide  a  deviation  from  a  best  fit 
straight  line  or  a  video  deviation  DEW.  After  focusing  (com- 


-72- 


pressing  in  range)  the  image  location  of  the  WC's  becomes  a  W 
in  the  azimuth  compression,  and  the  DEW  becomes  a  contributor 
to  range  smear. 

In  MSTRIP,  the  short  pulse  simulation,  an  input  RRMIN  is 
used  with  DEVR  to  give  range  smear  as 

smr  .  +  (rrmin)2 

(1  +  SLW  ) 

where  SLW  is  the  slope  of  the  W  vs.  G  fit.  Since  1  +  (SLW)^  is 
2 

the  esc  of  the  slope  angle,  dividing  by  it  provides  the  contri¬ 
bution  of  DEVR  in  the  direction  of  W  or  the  range  direction. 

6.4.2  Optimization  Choices 

The  second  history  in  Fig.  20  (A's)  is  a  good  illustration 

of  a  signal  history  which  can  be  optimized  with  respect  to  azimuth 
resolution.  Optimization  is  performed  in  the  following  way. 

First,  the  doppler  history  for  a  selected  target  (the  focus 
target)  is  screened  relative  to  the  minimum  (FSMIN)  and  maximum 
(FSMAX)  values  of  FS.  Dopplers  outside  this  range  are  not  used. 
Secondly,  each  possible  segment  that  contains  at  least  three  points 
and  consists  of  a  beginning  point,  an  ending  point,  and  all  points 
in  between,  is  checked  for  consistancy  in  sign.  Either  all  FS's 
are  positive  or  all  must  be  negative;  otherwise  this  segment  is 
not  used.  For  each  segment  that  remains,  the  azimuth  resolution, 
for  azimuth  optimization,  and,  if  RMS  optimization  is  desired, 
the  range  resolution,  for  the  focus  target  are  computed.  Table 
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IIF  presents  the  azimuth  resolute.',  table  for  the  target  with 
A-points  in  Fig.  20.  Data  for  this  figure  were  taken  from  FS 
and  Z  for  target  2  xn  Table  IIF.  The  smallest  value  in  this  table 
corresponds  to  the  best  possible  azimuth  resolution.  From  the 
table  we  obtain  9.678  as  the  minimum  azimuth  smear  (best  re¬ 
solution).  Tabular  output  is  organized,  so  that  each  row  cor¬ 
responds  to  a  fixed  number  of  points  in  each  segment.  The  first 
row,  one  entry,  corresponds  to  the  use  of  all  21  sample  points; 
the  second  row,  two  entries,  corresponds  to  segment  (1,  20)  and 
segment  (2,  21),  and  so  on  for  all  possible  segments.  (If  ***'s 

appear  in  an  entry,  the  corresponding  segment  was  excluded  by  the 

fch 

screening  process.)  Thus,  in  the  table,  an  entry  in  the  I 

th 

row,  from  the  top,  and  K  column,  from  the  left,  is  for  the 
segment  (Nl,  N2)  when 

N1  -  K  and  N2  -  21  -  I  +  !C 

Hi  nee  entry  9.678  is  for  segment  (15,21).  It  is  noted  that  a 
row  with  greater  than  13  entries  is  extended  around  onto  the 
next  line  and  an  extra  space  is  printed  after  an  extended  line. 

If  RMS  resolution  is  desired,  then  an  RMS  resolution  table 
is  generated.  RMS  optimization  is  computed  with  rules  similar 
to  those  for  azimuth  optimization.  The  only  difference  is  in  the 
table  used  for  choosing  the  best  resolution. 

After  optimization,  each  target  history  is  processed 
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(azimuth,  range,  and  RMS  resolutions  are  computed)  using  the  fre¬ 
quency  band  limits  and  focal  lengths  obtained  during  optimization. 

In  addition,  to  frequency  band  limits  (both  positive  or  both 
negative),  one  further  criteria  must  be  met  by  each  segment  of 
doppler  history.  That  is,  the  dopplers  are  required  to  be  mono- 
tonically  decreasing  with  increasing  Z.  Focal  length  at  the 
focus  target  is  directly  proportional  to  the  slope  cf  FS  vs.  Z. 

This  is  azimuth  focal  length.  It  is  assumed  for  subroutine  PRO¬ 
CESS  that  either  short  pulse  or  de-chirp  (range  compressed)  data 
is  used  as  an  input.  Targets  at  different  ranges  than  the  focus 
target  can  be  focused  only  at  a  longer  (for  greater  range)  or  shorter 
(for  near  range)  focal  length  that  is  proportional  to  range. 
Correcting  the  optimum  focal  length  by  the  ratio  of  range 
at  the  focus  target  to  the  range  of  the  target  being  processed 
produces  the  desired  change.  Doing  so  provides  the  essence  of 
a  tilted  plane  processor. 

If  chirp  pulses  (two  dimensional  signal  film)  are  to  be  pro¬ 
cessed,  then  subroutine  RNGPP.C  is  used  to  de-chirp  the  video 
signal  histories.  Subsequent  use  of  PROCESS  proceeds  as  discussed 
above.  Transparent  to  the  user  is  the  fact  that  RPMIN  is  computed 
in  using  RNGPRC  and  is  assumed  as  an  input  parameter  when  only 
PROCESS  is  used.  RRMIN  is  data  initialized  but  can  be  changed 
by  a  Type  II  parameter  input  to  any  desired  value. 
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6.4.3  Image  Location 

Images  are  located  in  azimuth  and  range.  A  focal  point 
for  image  location  is  obtained  by  linearly  fitting  W  vs.  FS 
and  G  vs.  FS.  W  is  range  location  of  signal  points  and  G  is 
azimuth  location  of  signal  points.  These  best  fit  straight  lines 
are  then  used  to  extrapolate  to  FS  ■  0.  FS  *  0  corresponds  to  the 
mean  focal  point.  Thus,  Wq  *  W(FS  *  0)  and  Gq  *  G(FS  *  0)  are 
the  linearly  predicted  coordinates  of  the  image  point. 

Two  important  distortions  in  image  location  are  due  to 
squint  angle  and  depression.  Since  the  nominal  squint  and  de¬ 
pression  are  known  a  priori,  geometrical  transformations  can 
be  performed  to  correct  for  these  distortions.  Corrections  for 
squint  can  be  accomplished  in  an  optical  processor.  Both  cor¬ 
rections  can  be  achieved  in  a  digital  system.  The  correction 
for  squint  is 

Gol  "  V  |cos  0cl  +  (W0  +  RMIN)  *  Sin  6c 

and 


W  .  *  W  *  cos  0 
ol  o  1  c1 

where  9c  is  the  squint  angle  to  center  of  beam,  Wq  and  Gq  are 
initial  location  of  focused  image,  and  adjust  from  range 

into  swath  to  total  range.  Gq  is  already  in  absolute  terms. 

The  effect  of  a  squint  correction  is  shown  in  Fig.  21  where 
A  is  the  input  image,  B  is  the  initial  output  image,  C  is  the 
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image  corrected  for  squint,  and  0  is  image  C  corrected  for  de¬ 
pression  angle.  The  purpose  of  the  first  correction  is  to  change 
B  by  erecting  the  image  so  that  equal  azimuth  inputs  tend  to  produce 
equal  azimuth  in  the  image.  The  second  correction  provides  a 
scale  change  in  the  range  direction  which  compensates  for  non- 
linearities  due  to  changes  in  depression.  This  latter  change 
is  effected  by 

U  »  G  , 
ol 

and 

V  -  W  ,/cos(DEPC) 

01 

where  DEPC  is  the  nominal  depression  angle  to  center  of  beam 
at  center  range  of  strip  swath. 

6.4.4  Geometrical  Distortion 

A  measure  of  geometrical  distortion  is  given  by  forming 
the  ratio  of  the  length  of  a  line  in  the  output  image,  to  length 
of  the  same  line  in  the  input  image.  For  example,  in  Fig.  6.17, 
Bj^/A^  is  a  measure  of  distortion.  Obviously  going  to  RSMOD 
C'1C2^1^2  improve  this  particular  ratio.  Similarly,  distor¬ 

tion  ratios  relative  to  the  D  image  (at  least  in  the  range  direc¬ 
tion)  will  be  better  than  for  C.  Since  target  inputs  to 
are  point  target,  it  is  a  simple  calculation  to  provide  every 
possible  ratio  for  the  distance  between  pairs  of  points  with 
respect  to  B,  C,  and  D  images. 
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6.4.5  Processor  Subroutines 


Imaging  in  the  RSMOD  simulation  is  accomplished  by  two  sub¬ 
routines.  To  compress  or  focus  video  history  a  range  focusing 
subroutine,  RNGPRC,  is  called  into  play.  This  subroutine  calls 
on  the  linear  fit  subroutines,  FIT  and  DEV,  to  obtain  range 
smear,  intermediate  range  and  azimuth  focussed  location,  and  to 
provide  a  contribution  of  video  non-linearity  to  azimuth  smear. 

Use  of  RNGPRC  is  implied  by  the  box  in  Fig.  10(b)  which  contains 
"Compress  Range  History  to  obtain  one  point  on  Azimuth  History". 
From  a  programming  sense  RNGPRC  fits,  most  naturally,  inside 
the  time  iteration  in  SIGGEN.  During  the  range  compression  for 
one  pulse,  video  spatial  frequencies  are  fitted  to  a  best  fit 
straight  line  with  respect  to  range.  Thus  range  smear  is  obtained. 
Dopplers  for  this  same  pulse  are  also  fitted  with  respect  to 
range,  to  obtain  a  smear  in  azimuth  due  to  non-linearity  of  doppler 
that  results  from  chirp  non-linearity  and  due  to  any  non-linear 
sampling  rates  in  range.  Similarly  the  focused  point  in  range, 
where  video  spatial  frequency  goes  to  zero,  is  also  obtained. 

From  the  linear  fit  of  doppler  vs.  range,  the  mean  doppler  to 
be  associated  with  this  range  (at  zero  video)  io  easily  computed. 

After  signal  generation  and/or  range  processing,  the  sub¬ 
routine  PROCES  is  called  to  obtain  final  images  and  smear  as 
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described  earlier.  Figure  23  is  a  generalized  flow  diagram 
for  this  latter  subroutine. 

6.5  Output 

There  are  numerous  integer  and  logical  switches  that  con¬ 
trol  output.  For  the  purposes  of  this  discussion,  we  illustrate  a 
useful  selection  of  outputs  as  shown  in  Table  II  and  III.  The 
order  of  outputs  are  fixed. 

Considering  individual  parts  of  the  output  we  have  first 
the  run  number,  RUN  NO.  1  in  this  case,  followed  by  a  line  headed 
with  &  PARAMS.  This  latter  line  indicates  parameter  values 
selected  as  TYPE  II  inputs  for  this  run.  Following  the  params 
line  is  a  table  of  roll,  pitch,  yaw,  xs,  ys,  zs,  vxs,  vys,  and  vzs 
which  are  target  frame  motions  provided  through  the  use  of  the 
subroutine  SHIPMO.  The  table  is  100  items  long  and  corresponds 
to  a  10  second  period  of  logged  ship  motion.  Underneath  this  table 
is  a  list  of  parameters  used  in  developing  the  tabular  entrees. 

As  labeled  this  motion  history  came  from  Buckley  logs, 
coded  14A1.  Also,  the  center  item,  item  number  50,  corresponds 
to  a  GMT  (Greenwich  Mean  Time)  of  14  hrs.,  14  min  and  53.0  seconds. 

Table  IIC  continues  the  example  output  with  three  lines 
about  doppler  history.  For  each  of  the  three  targets,  the 
number  of  useable  parts  (segments)  and  the  largest  part,  with 
the  beginning  point  and  ending  point,  are  given.  Following  these 
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Figure  22.  PROCES  Generalized  Flow  Diagram. 
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SUMM*uy  OF  P*PA«£T£P  CHANGES: 


results  of  azimuth  processing  are  t*?o  items,  FMBEAM  and  FMRATE, 
for  each  target  which  are  FS  at  mid  beam  and  linear  fit  slope 
of  FS  vs  2.  These  particular  values  are  computed  for  no  motion 
situations  in  which  the  aircraft  flies  perfectly.  FMBEAM  and 
FMRATE  are  then  used  by  the  impulse  response  calculations  to  form 
a  standard  filter  to  compress  the  actual  doppler  history  observed 
under  the  prescribed  mocion  of  a  simulation. 

A  complete  listing  of  Type  II  inputs  is  shown  at  the  top 
of  Table  HD*  In  this  particular  listing,  the  value  -  0.69902160E-76 
corresponds  to  a  core  constant,  computer  dependent,  and  a  storage 

location  which  is  not  being  used.  For  this  example  table  only, 

three  targets  are  being  used  but  fifteen  locations  are  reserved 
for  each  target  oriented  parameter. 

Then  a  complete  history  of  dopplers  for  every  target  along 
with  a  number  of  interesting  parameters  is  printed.  The  following 
definition  hold: 

T  -  time  from  beginning  of  run 

VX  -  aircraft  velocity  x  direction 

VY  -  aircraft  velocity  Y  direction 

VZ  -  aircraft  velocity  Z  direction 

R  -  instantaneous  slant  range  between  target  point 

and  phase  center  of  radar 
SINTHA  -  sin  of  actual  squint  angle  to  target 

SINDEP  -  sin  of  depression  angle  to  target 

RDOT  -  Range  rate  or  beam  velocity  (positive  is  closing) 

PSRATE  -  Phase  shifter  rate  in  ft/sec 

FTEMP  -  Temporal  frequency  of  doppler  in  cycles/sec  after 
phase  shifting  and  clutterlock  correction 
FILMSP  -  film  speed  in  feet/sec. 

FSPAT  -  spatial  frequency  in  cycles/ft. 

G  -  Azimuth  position  in  feet 

W  -  Range  position  in  feet 

Z  -  Position  along  doppler  trace  in  feet 


Due  to  target  motion  we  note  that  RDOT  and  hence  FSPAT  for  target 
1  changes  sign  and  FSPAT  is  not  monotonically  decreasing.  Also 
FSPAT  for  target  1  is  much  lower  than  for  target  2,  the  focus 
target.  As  a  consequence  target  number  1  has  no  useable  parts 
as  indicated  in  Table  IIC. 

Following  the  doppler  history  print  out,  in  Table  1IF,  we 
have  the  azimuth  resolution  table  for  the  focus  target  and  a 
summary  of  the  best  fit  straight  line  which  results  from  azimuth 
focusing.  The  range  and  rms  resolutions  which  correspond,  along 
with  average  values,  are  also  given.  Then  a  table  of  target  posi¬ 
tions  G  and  W  (direct  from  signal  films),  G1  and  G2  (compensated 
for  squint  -  intermediate) ,  and  U  and  V  (compensated  for  depres¬ 
sion  -  final  film)  are  presented.  Distortion  ratios  for  each  of 
these  three  output  locations  are  then  shown  on  the  next  page  of 
print  out. 

For  a  quick  look  and  for  perspective  orientation  in  smear 
and  final  location,  the  G  and  W  positions  (or  alternatively  the 
G1  and  Wl)  positions  can  be  plotted  as  shown  in  Table  IIH.  A 
scale  in  feet  per  row  and  feet  per  column  (character  position  in 
a  row)  and  a  heading  of  run  number  and  target  information  heads 
this  print  out.  Numerals  are  then  used  to  plot  input  target 
location  with  letters,  correspondences  shown  by  the  list  on  the 
right  side,  to  show  image  locations.  Multiple  letters  like  BB 
indicate  how  much  azimuth  smear  has  occurred. 
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Finally,  the  last  page  of  print  out  contains  a  summary 
by  run  number  of  the  specific  parameters.  Type  II,  used  in  each 
run. 

Table  HI  shows  an  alternative  in  output.  In  this  parti¬ 
cular  case,  amplitude  of  target  return  and  the  roll  rate  history 
of  the. target  were  selected  for  output.  Amplitude  is  in  watts. 

No  target  motion  was  allowed  for  this  table  and  as  a  consequence 
the  entire  synthetic  aperture  for  each  target  makes  up  a  useable 
part.  For  the  param  lii.e,  Table  4A,  LNEWP  =  F  will  stop  the 
program  with  $END,  LSIG  =  T  will  cause  a  new  signal  film  to  be 
generated,  and  LAMPR  =  F  signifies  that  no  new  target  patterns 
are  to  be  read  in.  Actually  LAMPR  was  set  true  with  Type  I  input, 
then  it  is  automatically  set  false  before  Type  II  inputs  begin. 
Hence,  LAMPR  *  F  in  the  params  line  is  redundant. 

6.6  Conclusions 

With  respect  to  the  RSMOD  model,  all  of  the  goals  set  out 
in  SAR  DIGITAL  Simulation  Program  Plan  have  been  reached.  Com¬ 
pilation  of  the  entire  computer  model  is  essentially  complete. 

All  but  the  final  checking  and  linking  together  is  done.  With 
the  exception  of  chirp,  chirp  processing  and  the  MTF  functions 
all  RSMOD  functions  have  been  used  to  obtain  inputs  for  the 
impulse  response  calculations. 

There  are  two  omissions  in  the  presently  coded  program. 

First  a  function  and  a  call  statement  are  required  to  realize 
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the  effect  of  a  video  amplifier.  This  omission  can  be  easily 
fixed  since  the  F.F-Amplifier ,  which  is  coded,  is  almost  identical 
to  the  video  amplifier  representation.  Secondly,  the  signal 
recorder  in  the  present  version  does  not  contain  a  transfer 
function  to  attenuate  signal  strengths  according  to  the  input 
of  a  characteristic  function.  For  example,  a  characteristic 
function  is  required  in  an  optical  recorder  representation  whenever 
a  signal  level  (volts)  is  transformed  into  light  intensity  on 
the  face  of  a  CRT  or  when  this  intensity  is  transformed  by  a 
film  characteristic  into  a  gray  scale  on  the  film.  A  mathe¬ 
matical  representation  of  each  of  these  functions,  amplifier  or 
characteristic  function,  is  available  but  neither  have  been 
coded  into  the  program.  Each  of  the  characteristic  functions 
can  be  input  as  a  table  obtained  from  calibration  experiments 
or  from  equipment  specifications. 
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7,  IMPULSE  RESPONSE  GENERATION 


7 . 1  Introduction 

A  major  item  within  the  scope  of  the  SAR  Imaje  Quality  Analysis  program 
is  the  development  of  a  digital  computer  model  of  a  SAR  radar  system  which 
is  to  be  utilized  for  determining  the  effects  of  radar  parameter  variation  on 
image  quality.  The  model  has  been  constructed  so  as  to  determine  the 
response  of  an  airborne  SAR  system  in  the  presence  of  realistic  errors 
(e.g. ,  imperfect  motion  compensation,  transmitter  chirp  nonlinearities) 
under  actual  aerodynamic  flight  conditions  characteristic  of  high  performance 
aircraft.  An  assessment  of  image  quality  is  to  be  obtained  in  part  by  an 
analysis  of  the  response  of  this  model  to  point  target  inputs  (i.e.,  the  impulse 
response) .  This  section  will  describe  the  software  module  which  has  been 
developed  to  accept  SAR  simulation  outputs  and  produce  an  actual  point 
target  image.  At  present,  the  model  is  one-dimensional  and  is  capable 
of  processing  either  in  azimuth  or  range.  Extensions  to  two-dimensionai 
processing  is  straightforward  and  could  be  implemented. 

7. 2  Theory  of  Operation 

The  SAR  model  is  designed  to  be  a  "high-fidelity"  digital  simulation  of 
a  side-looking,  chirp-pulsed,  synthetic  aperture  radar.  In  either  the  range 
or  azimuth  direction,  and  under  the  usual  far-field  assumptions,  the  signal 
pulse  returned  from  a  point  target  is  essentially  of  the  form 
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where, 


r(x)  =  a(x-xQ)  cos  (4>(x-x^)) 


<7-l) 


0(x)  =  u  X 
0 


ax 


V-2) 


and  a  is  a  parameter  proportional  to  the  linear  FM  rate  of  the  pulse.  For 
azimuth  processing  a  is  proportional  to  l/AR  where  A  is  the  radar  wave¬ 
length  and  R  is  the  point  target  range.  For  range  processing,  a/2  is 
equivalent  to  the  FM  rate  of  the  chirp  pulse.  The  term,  is  the  mid¬ 
band  frequency  of  the  chirped  pulse  in  the  range  processing  case,  or  is 
the  doppler  offset  caused  by  antenna  squint  in  the  azimuth  case.  The 
processing  required  in  both  cases  is  identical  so  that  the  software  described 
in  this  section  is  applicable  to  processing  of  either  azimuth  or  range  data. 

The  variable  x  in  the  above  equations  can  refer  to  time,  slam-range  distance, 
or  film  spatial  coordinates,  depending  upon  the  scale  factor  inherent  in  x. 
Presently  the  variable  used  in  the  impulse  response  programs  refers  to  slant 
range  distance.  The  constant,  xq,  represents  the  target  position.  The 
amplitude  term,  a(x),  is  a  function  of  antenna  pattern  and  target  cross- 
section  and  presently  is  assumed  to  be  a  real  quantity. 

The  processing  of  pulses  of  the  form  (7-1)  into  point  target  images  can  be 
accomplished  by  a  matched  filtering  operation.  That  is,  the  signal  in  (7-1)  is 
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input  to  a  linear,  time-invariant  system  whose  impulse  response  is  given 
by  h(x)  =  r(-x),  where  to  bar  (")  indicates  a  complex  conjugate  operation. 
The  resulting  image  is  represented  as 


g(x)  =  r(x)*h(x)  =  /  r(r)h(x-T)dr. 


The  operation  given  in  (7-3)  is  most  often  performed  in  an  optical  analog 
fashion,  but  may  be  performed  digitally  as  in  the  present  model.  There 
is  generally  a  finite  aperture  over  which  the  convolution  can  take  place, 
such  that  g(x)  is  more  correctly  given  as 


/oo 
r< 


'X-r  +  x 


r(r)h(x-r)  rect 


Under  realistic  operating  conditions,  system  errors  introduce  undesirable 
perturbations  into  the  ideal  received  signal  as  given  in  (7-1).  Thus,  instead 
of  processing  r(x)  one  must  work  with  a  signal  represented  as 


r'(x)  =  a’(x-xQ)  cos  (^'(x-x^)  (7-5) 

which  is  at  best  only  an  approximation  to  the  signal  given  in  (7-1).  It  is  an 
analysis  of  the  resultant  signal, 
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g’(x)  =  r'(x)*h(x) 


(7-6) 


that  will  produce  an  assessment  of  image  quality.  To  generate  g'(x)  one 
musf  therefore:  1)  simulate  the  SAR  system  in  order  to  generate  arti¬ 
ficially  r'(x);  2)  generate  analytically  an  appropriate  h(x);  and  3)  perform 
the  convolution  operation.  Each  of  these  tasks  will  be  discussed  in  detail 
in  relation  to  the  current  SAR  simulation  program. 

7.2. 1  Received  Signal  Generation.  The  SAR  simulation  package 
(RSMOD)  is  structured  so  as  to  provide  a  set  of  equal  time  samples  of 
a'(x)  and  of  0'(x),  where  <j>'(x)  represents  the  instantaneous  frequency  of 
the  returned  pulse.  Due  to  the  possible  presence  of  motion  errors  there 
may  no  longer  be  a  one -one  correspondence  between  time  and  spatial 
coordinates.  Therefore  a  spatial  coordinate  v?  ue  is  supplied  with  each 
equal-time  sample  of  signal  amplitude  and  frequency.  Knowledge  of 
the  sample  values  and  their  spatial  positions  allows  one  to  apply  inter¬ 
polation  techniques  to  obtain  equidistant  amplitude  -  frequency  samples. 

Only  equidistant  samples  can  be  efficiently  processed  in  the  digital  environment. 
In  addition,  the  interpolation  step  gives  one  the  opportunity  to  generate  a 
sufficient  number  of  samples  to  satisfy  the  Nyquist  sampling  criterion  for 
the  r'(x)  pulse  to  be  generated  (#  samples  >  2-TBP  or  the  pulse,  where 
TBP  equals  the  duration  of  the  pulse  times  its  bandwidth.)  The 
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simulation  program  can  supply  far  fewer  samples  of  a'(x)  and  i'(x)  than 
are  required  to  represent  r'(x/,  since  generally  a'(x)  and  *'(x)  vary  slowly 
with  respect  to  x  compared  to  r'(x),  which  can  be  a  relatively  high  frequency 
pulse. 

A  Chebyshev  interpolation  scheme  has  been  chosen  for  this  application. 

Let  f(x  ),  f(x  )  ...,  f(x  )  be  samples  (not  necessarily  equidistant)  of  f(x), 
i  £4  n 

where  f(x)  could  represent  either  a’(x)  or  *'(x).  The  problem  is  to  determiie 
the  coefficients  of  the  polynomial 

n  1 

P(x>=£  ex'  ,  (7-7) 

i=l 

such  that 

"  2 

/  ,  (f(x.)  -  p(x.))  =  minimum.  (7-8) 

i=l  1  1 

Once  the  coefficients  c^,  i=l,  ...  n  are  found  one  can  produce  as  many 
equidistant  samples  of  p(x)  as  desired  simply  be  evaluating  the  polynomial. 

A  Chebyshev  interpolation  involves  expansion  of  the  polynomial  p(x)  in  terms 
of  its  Chebyshev  expansion.  The  Chebyshev  expansion  of  p(x)  gives  a  much 
better  indication  of  the  accuracy  of  the  curve  fit  than  the  coefficient  vector 
of  the  polynomial  itself.  If  the  specified  degree  of  the  polynomial  is  too 
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high,  the  last  terms  cf  the  Chebyshev  expansion  will  be  uniformly  small 
compared  to  lower  degree  coefficients.  The  degree  may  be  reduced  by 
the  number  of  small  trailing  coefficients  without  unduly  enlarging  the 
overall  error.  This  possibility  is  checked  for  automatically  by  the 
present  interpolation  program.  The  reader  is  referred  to  the  IBM  System/360 
Scientific  Subroutine  Package  Programmer's  Manual  for  further  details  con¬ 
cerning  the  interpolation  scheme.  The  subroutines  APCH  and  APFS,  are 
used  by  the  subroutine  SFIT  to  gererate  the  Chebyshev  coefficients.  The 
subroutine  named  CNPS  is  used  to  evaluate  the  pohmomial  p(x)  to  obtain 
the  new  function  values . 

Following  the  interpolation  procedure  one  has  M  equidistant  samples 
of  a'(x)  and  o'(x)  with  which  to  generate  M  samples  of 

r'(x)  =  a'(x)  cos  (7-9) 


That  is,  the  values 


and 


a'(knx),  k  =  1,  . . . ,  M 


i'(kox),  k  =  1,  ...,  M 
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are  stored  in  computer  memory  for  use  in  generating  the  received  signal, 
r'(x).  Obviously  one  requires  a  method  of  generating  samples  of  <*>'(x) 
given  samples  of  ^'(x).  Recall  the  fundamental  theorem  of  calculus.  If 
f(x)  is  continuous  for  a  <  x  <  b,  and  if  f  (x)  is  any  anti-derivativo  of  f(x) 
on  this  interval  (i.e. ,  F  (x)  =  f(x)  for  a  <  x  <  b),  then 


(7-10) 


We  can  assume  that  *'(x)  is  continuous  over  the  interval  of  interest  and 
that  <f>'(x)  is  its  anti-derivative  on  interval  (by  definition).  Therefore 
by  the  Fundamental  Theorem, 


’(y)dy  =  *'(x)  -  *’(0) 


(7-11) 


or  rearranging, 


*'(x)  =  <f>’(0)  +  j’*'(y)dy.  (7-12) 


In  terms  of  digital  quantities, 

k 

*'(kAx)  -  A  +^i'(jnx),  K=  1,  ....  M,  (7-13) 

j=l 
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».  '  **e  A  =  $'(0),  or  the  phase  of  the  signal  immediately  upon  the  target 
entering  the  beam.  This  constant  is  chosen  such  that  the  phase  of  the 
signal  is  zero  at  midbeam,  and  the  reference  signal  (described  later) 
is  derived  accordingly.  Having  produced  <$ (k Ax) ,  k  =  1,  . . .,  M  and  using 
a'(kAx)  from  the  interpolation  procedure,  the  received  signal  samples 
are  generated  as 

r'(kAx)  =  a*(kAx)  cos  ($'(kAx)),  k=l,  ....  M.  (7-14) 

The  values  of  r'(kAx)  are  equivalent  to  samples  taken  at  equidistant  inter¬ 
vals  (and  satisfying  the  Nyquist  criterion)  in  the  film  or  signal  plane. 

7,2.2  Reference  Function  Generation.  The  reference  function  refers 
to  the  impulse  response  of  the  linear  time-invariant  system  used  to 
compress  the  received  pulses.  In  the  present  case  it  is  of  the  form 

h(x)  =  exp{j(«ox  +ax2)| .  (7-15) 

A  complex  reference  function  is  used  so  that  the  carrier  («q)  inheT  mt 
in  the  real  received  si^al  can  be  easily  removed  and  to  overcome  problems 
caused  by  conjugate  images.  In  the  computer,  h(x)  (or  rather  h(kAx))  may 
be  generated  analytically  with  knowledge  of  and  a.  Several  options  have 
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been  programmed  as  method?  ■  I  choosing  uQ  and  a.  They  are:  1)  Using 
the  theoretically  "correct"  and  a  for  each  target.  These  values  are 
entered  into  the  subroutine  as  arguments  and  can  be  supplied  by  the 
user,  or  alternatively,  calculated  by  the  SAR  simulation  program  itself. 

2)  Resorting  to  an  autofocus  mode  whereby  and  a  are  estimated  using 
the  received  doppler  frequencies.  3)  Using  a  previously  generated 
reference  function.  In  this  case,  the  reference  function  used  for  com¬ 
pressing  the  previous  target  is  just  re-used,  thereby  reducing  the 
amount  of  processing  time  required.  A  user  must  be  cautious  if  relying 
upon  this  third  option  in  ensuring  that  the  reference  function  is  the  correct 
one.  For  instance,  if  one  is  azimuth  processing,  then  only  a  string  of 
targets  of  constant  range  can  be  processed  with  the  same  reference 
function  since  a  is  a  function  of  range  for  azimuth  processing. 

Processing  option  number  two  above  requires  additional  explanation. 

In  this  case  an  estimate  must  be  made  of  wq  and  a  using  the  received  samples, 

i'(kAx).  The  estimated  midbeam  frequency  «  is  simply  set  equal  to  the 

o 

time  sample  of  $'(x)  which  comes  nearest  to  the  midbeam  position.  An 
estimate  of  o,  denoted  «,  is  generated  using  the  1st  order  Chebyshev 
coefficient  derived  in  the  interpolation  process.  The  resultant  FM  rate, 


-,  is  equivalent  to  a  least-squares  estimate  of  the  FM  rate  of  the  incoming 
2 

pulse.  In  this  manner  the  system  tends  to  optimize  the  match  between 
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r'(x)  and  h(x),  thus  optimizing  the  focus  for  each  target  it  processes.  Since 
the  estimates  a»o  and  a  are  already  by-products  of  the  interpolation  process 
the  autofocus  option  requires  no  additional  computation  time  over  that  required 
for  the  theoretically  "correct”  focus  option. 

7.2.3  Performing  the  Convolution.  The  convolution  required  to  compress 
the  SAR  signals  must  be  performed  digitally.  The  form  of  the  operation 
originally  given  by  (7-3)  is  thus  of  the  form 

qo 

g(kAx)  =  ^  r'(iAx)h((k-i)Ax),  (7-16) 

i=— 00 

where  it  is  assumed  that  r’(kAx)  and  h(kAx)  are  both  zeio  for  k<  1  and 
K  >  M.  The  number  of  computations  involved  in  performing  the  di^'^ete 
convolution  of  (7-16)  can  be  significantly  reduced  by  use  of  the  Fast 
Fourier  Transform  (FFT)  algorithm.  The  steps  required  in  this  process 
are  listed  as  follows: 

1)  Pad  at  least  M  zeros  onto  the  ends  of  the  arrays  containing 
v'(kAx)  and  h(kAx)  to  insure  that  an  aperiodic  convolution 
will  take  place. 

2)  Take  the  Discrete  Fourier  Transform  (DFT)  of  the  padded  sig¬ 
nals  r'(kAx)  and  h(kAx)  using  the  FFT  algorithm  to  save  com¬ 
putation  steps.  Denote  the  transformed  data  as  R'(kAf)  and 
H(kAf),  respectively. 
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3)  Multiply  R'(kAf)  by  H(kAf)  on  a  point-by-point  basis, 
k=l,  . . M. 

4)  Compute  the  inverse  DFT  of  the  product  R'(kAf)-H(kAf). 

This  yields  the  desired  convolution  output. 

Step  1)  above  is  required  since  multiplication  of  DFTs  corresponds 
to  convo  '.ion  of  periodic  functions,  not  aperiodic  ones.  A  convolution 
of  pej  .die  functions  can  also  be  regarded  as  a  circular  convolution 
in  which  values  of  the  kernal  that  are  shifted  from  one  end  of  a  period 
are  circulated  into  the  other  end,  thus  introducing  what  has  been  termed 
interperiod  interference.  By  the  proper  padding  with  zeros  a  periodic 
convolution  can  be  used  to  compute  an  aperiodic  convolution.  This 
is  done  by  adding  to  each  cycle  of  both  functions  a  string  of  zero  values 
of  sufficient  width  to  provide  a  kind  of  "grace  band"  during  which  the 
interperiod  interference  assumes  zero  numerical  value. 

In  performing  step  2)  some  computation  time  can  be  saved  due  to  the 
fact  that  r'(kAx)  is  a  real  valued  function.  Transforming  of  a  real  array 
of  length  2M  can  be  accomplished  by  transforming  a  complex  array  of 
length  M  and  then  rearranging  the  resulting  array.  The  point-by-point 
multiplication  of  step  3)  is  then  accomplished  by  multiplying  the  positive 
frequency  terms  of  H(kAf)  by  the  terms  of  R'(kAf)  and  setting  the  remaining 
terms  of  the  product  equal  to  zero.  This  procedure  follows  from  tacitly 
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assuming  that  the  spectrum  of  h(kAx)  contains  no  negative  frequency 
energy.  Thus  the  present  algorithm  assumes  the  sufficient  azimuth  offset 
case  whereby  non-positive  frequencies  are  absent.  This  procedure  can  be 
easily  modified  to  account  for  negative  frequencies.  The  presence 
of  the  non-positive  frequency  spectrum  would  create  conjugate  images 
and  D.  C.  term  interference  in  the  current  processing  operation.  In 
any  case,  step  4)  is  a  straightforward  complex  transform  operation. 

Ideal  system  performance  occurs  when  the  reference  function  is 
perfectly  matched  to  the  incoming  pulse,  as  is  the  case  when  h(x)  =  r(-x). 

In  this  case  H(kAf)  =  R(kAf),  the  conjugate  of  R(kAf).  Therefore,  for 
a  given  reference  function,  the  ideal  system  response  can  be  obtained  by 
multiplying  H(kAf)  by  its  conjugate  and  then  inverse  transforming  the 
result.  This  procedure  is  included  in  the  present  system  as  a  user 
option.  If  the  user  specifies  the  use  of  this  option,  the  program  derives 
the  ideal  system  response  for  the  reference  function  of  every  target  response 
generated. 

The  data  generated  by  the  point  target  response  software  is  utilized 
primarily  in  two  ways;  visual  display  and  digital  analysis.  For  a  visual 
assessment  Calcomp  and  computer  printer  plots  are  produced  of  signal 
frequency,  signal  phase,  signal  recording  (zone  plate  section),  reference 
phase,  and  target  response.  Target  response  may  be  displayed  as  an 
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amplitude  plot,  a  logarithm  of  intensity  plot,  or  an  expanded  log-intensity 
plot.  Ideal  target  response  can  be  generated  and  displayed  in  an  identical 
fashion.  Digital  analysis  of  the  point  target  response  is  to  be  performed 
via  the  Modulation  Transfer  Function  (MTF).  To  this  end  the  MTF  is 
generated  and  returned  to  the  calling  program  in  an  auxiliary  array. 

The  MTF  is  the  modulus  of  the  Fourier  transform  of  the  point  target 
intensity  response.  After  multiplication  of  R'(kAf)  and  H(kAf)  one 
has  DFT  of  amplitude  response.  Therefore,  one  must  take  an  inverse 
FFT  to  get  amplitude  response,  square  the  resultant  data,  and  then  perform 
one  additional  FFT. 

All  of  the  processing  and  display  options  mentioned  above  are  con¬ 
trolled  by  the  user  (calling  program)  through  the  control  parameters  in 
the  argument  list.  The  impulse  response  software  package  (IRMOD)  is 
in  the  form  of  a  subroutine  which  may  be  called  by  the  simulation  package 
(RSMOD).  In  addition,  through  the  use  of  special  purpose  interface  programs, 
the  impulse  response  package  can  be  operated  offline.  The  next  section  will 
describe  the  software  in  detail. 

7.3  Target  Response  Software  Description 

Table  IV  lists  and  identifies  the  arguments  passed  to  and/or  received 
from  subroutine  PROCES,  the  primary  subroutine  making  up  IRMOD. 


Appendix  A  contains  a  detailed  flow  diagram  of  subroutine  PROCES. 
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TABLE  IV 

SUBROUTINE  PROCES  ARGUMENT  LIST 


ARGUMENT 

FUNCTIONAL  DESCRIPTION 

XPOINT 

A  real  *4  array  containing  the  spatial  positions  of 
ihe  doppler  history  samples. 

FREQU 

A  real  *4  array  containing  the  doppler  frequency 
of  each  of  the  doppler  history  samples. 

AMPL 

A  real *4  array  containing  the  amplitude  of  each 
of  the  doppler  history  samples. 

NPTS 

An  integer *4  location  containing  the  number  of 
doppler  history  samples  and  hence  the  number 
of  elements  in  each  of  the  three  doppler  sample 
arrays.  NPTS  must  be  <21  at  present. 

IPOW 

An  integer *4  location  containing  the  power  with 
which  to  attempt  the  Chebyshev  curve  fitting 
procedure.  IPOW  must  be  less  than  or  equal 
to  10. 

NPEXP 

An  integer*4  location  containing  an  exponent  of 
two  and  describing  the  number  of  samples  to 
be  generated  within  the  digital  processing 
aperture.  The  number  of  samples  is  equal  to 
2**nPEXP. 

IOCNTL 

An  integer*4  location  containing  all  option  speci¬ 
fications  in  coded  form.  Currently  the  13  least 
significant  bits  represent  I/O  and/or  processing 
options.  A  one  in  a  given  bit  location  indicates 
that  the  associated  option  is  to  be  in  effect.  The 
following  is  a  list  of  control  bit  assignments,  be¬ 
ginning  with  the  least  significant  bit  and  progress] 
toward  the  most  significant  bit: 
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ARGUMENT 


TABLE  IV  (Continued) 


FUNCTIONAL  DESCRIPTION 


(LSB) 


l)  set  for  obtaining  a  listing  of  all  original 
doppler  samples,  and  all  Chebyshev 
interpolation  coefficients. 


2)  set  if  all  computer  printer  plots  are 
to  be  Calcomp  plotted  also 


3)  set  to  obtain  a  plot  of  interpolated 
signal  amplitude 


4)  set  to  obtain  a  plot  of  interpolated 
signal  frequencies,  (doppler  history) 


5)  set  to  obtain  a  plot  of  the  signal  phase 
(the  integral  of  frequency) 


6)  set  to  obtain  a  plot  of  the  signal  itself 
(zone  plate  section) 


7)  set  to  obtain  a  plot  of  the  reference 
signal  phase. 


8)  set  to  obtain  a  plot  of  the  point  target 
amplitude  response 


9)  set  to  obtain  a  plot  of  the  logarithm  of 
target  intensity  response,  (truncated 
below  -40  db  and  normalized  to  0  db) 


10)  set  to  obtain  an  expanded  plot  of  the 
return  above  -40  db.  (Log  intensity) 
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TABLE  IV  (Continued) 


ARGUMENT  FUNCTIONAL  DESCRIPTION 


11)  set  to  compute  and  obtain  a  plot  of  the 
ideal  target  amplitude  response. 

12)  set  to  obtain  a  plot  of  the  logarithm 
of  the  ideal  target  intensity  response. 

13)  set  to  obtain  an  expanded  plot  of  the 
ideal  target  log-intensity  response. 

14)  unused 


LFOCUS  =3  implies  that  the  reference  function 
generated  for  the  previously  processed 
target  is  to  be  used. 

FMBEAM  A  real *4  location  containing  the  user  specified 

midbeam  frequency  to  be  used  in  reference 
function  generation  if  the  theoretically  "correct" 
focus  option  is  to  be  used.  (i.e. ,  if  LFOCUS=l) 

FMRATE  A  real*4  location  containing  the  user  specified 

linear  FM  rate  to  be  used  in  reference  function 
generation  if  the  theoretically  "correct"  focus 
option  is  to  be  used,  (i.e.,  if  LFOCUS=l) 

RESP  A  real *4  array  area  which,  upon  return ,  will 

contain  the  MTF  of  the  point  target  response. 
The  array  must  be  of  length  2**(NPEXP+1). 
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TABLE  IV  (Continued) 


ARGUMENT 


FUNCTIONAL  DESCRIPTION 


IER 


An  integer*4  location  which,  upon  return,  will 
contain  an  error  indication,  (i.e.,  iflER=2  the 
compression  of  the  signal  history  v*  .is  unsuccess¬ 
ful.) 


S-iy-.H.1,  »v.  " 


-:•=  ■  -u'  --  •  '  - 


Allowing  the  reader  to  use  Table  IV  and  the  Appendix  as  guides 
this  section  will  give  an  in-depth  description  of  subroutine  PROCES, 
the  primary  routine  in  IRMOD.  The  remaining  subroutines  are  general 
purpose  routines  that  are  largely  self-explanatory  and  will  not  be  described 
further  herein. 

Upon  entry  into  PROVES  a  set  of  program  constants  are  initialized  to 
correspond  with  the  user  specified  digital  aperture  size  (via  NPEXP) . 

Then  the  program  checks  to  insure  that  the  number  of  doppler  history 
samples  (NPTS)  does  not  exceed  21,  and  that  the  specified  power  of  the 
Chebyshev  curve  fit  does  not  exceed  10.  If  either  of  these  conditions  occur 
an  appropriate  error  message  is  printed  and  control  is  returned  to  the 
calling  program  with  the  error  indicator  set  (IER--2) .  The  next  step  is  to 
decode  the  I/O  control  parameter.  This  is  done  by  extracting  each  of  the 
13  bits  used  in  IOCNTL  and  storing  each  bit  for  later  reference  in  the  13 
elements  of  the  array  IPLT,  Next  the  Chebyshev  curve  fit  power  is  stored 
temporarily  in  locations  LPOW  and  MPOW.  Thus,  the  value  of  the  curve 
fit  power  for  either  the  amplitude  or  frequency  interpolation  procedure 
can  be  modified  later  if  the  accuracy  of  che  fit  warrants  it,  without 
modifying,  the  location  (IPOW)  passed  to  the  routine  by  the  user. 
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The  program  next  performs  the  Chebyshev  curve  fitting  procedure 
on  the  doppler  frequency  samples  supplied  by  the  calling  program.  The 
curve  fit  is  done  by  subroutine  SFIT,  which  returns  the  actual  power 
of  the  fit,  the  Chebyshev  coefficients,  data  scale  and  shift  factors,  and 
an  error  indicator.  The  scale  and  shift  parameters  (XD  and  XO)  are  used 
to  modify  the  doppler  sample  positional  data  (XPOINT)  such  that  the  curve 
fit  takes  place  on  the  new  interpolated  values.  If  the  curve  fit  was 
unsuccessful  (e.g. ,  the  required  minimum  error  was  not  obtained)  the 
error  parameter  will  have  the  value  "2"  upon  return  and  the  subroutine 
PROCES  will  print  an  appropriate  error  message  and  return  to  the 
calling  program. 

Before  performing  a  curve  fit  to  the  amplitude  data  the  program 
checks  to  see  if  the  amplitudes  associated  with  the  doppler  history  samples 
are  uniform.  If  the  amplitudes  are  all  equal  then  the  completed  cuv.e 
fitting  procedure  is  bypassed  and  a  logical  variable  named  AMPCHK  is 
set  to  true.  Otherwise,  a  call  to  SFIT  is  made  as  in  the  frequency  fit 
case  and  AMCHK  is  set  to  false.  The  power  of  the  amplitude  fit  and 
the  associated  coefficients  are  stored  in  separate  locations  than  those 
used  for  the  frequency  curve  fit.  If  the  amplitude  samples  are  all 
equal,  their  value  is  retained  in  location  SAVE  for  later  reference. 


Next,  the  linear  FM  rate  of  the  reference  function  is  defaulted  to 
that  predicted  by  the  1st  order  Chebyshev  (frequency  fit)  coefficient. 

Then,  if  LFOCUS=l  (theoretical  focus)  has  been  specified,  the  linear 
FM  rate  is  reset  to  that  value  specified  by  the  user  via  the  argument 
FMRATE  (see  Table  IV).  (Note  that  the  variable  ALPHA  represents 
one-half  the  FM  rate.)  The  same  procedure  is  followed  for  the 
reference  function  midbeam  frequency,  FMB.  The  midbeam  frequency 
is  first  defaulted  to  the  frequency  of  the  time  sample  which  occurs  nearest 
the  midbeam  position,  and,  if  LFOCUS=l  is  reset  to  that  value  specified 
by  the  user  via  the  subroutine  argument  list.  Thus,  the  two  parameters 
which  specify  the  reference  function,  the  linear  FM  rate  and  the  midbeam 
frequency,  are  assigned  according  to  the  user  specified  focus  option. 

The  program  next  enters  a  do-loop  in  which  the  interpolated  equally 
spaced  samples  of  frequency  and  amplitude  are  evaluated,  the  frequency 
is  integrated  to  produce  instantaneous  phase,  and  the  reference  function 
values  are  generated  —  all  within  the  same  loop.  In  addition,  the  values  of 
the  equally  spaced  sample  positions  are  generated  by  adding  the  sample 
spacing  to  the  previously  generated  sample  position.  Once  the  new  sample 
position  is  known  each  time  through  the  loop,  that  position  is  converted  to 
the  range  -1,  +1  using  the  scale  and  shift  factors  previously  returned 
from  SFIT  (i.e.,  XD  and  XO).  Using  this  converted  positional  value, 


the  frequency  and,  if  required,  the  amplitude  of  each  new  sample  are 
generated  by  calls  to  subroutine  CNPS.  (If  logical  variable  AMPCHK= 
.TRUE. ,  then  the  call  to  CNPS  can  be  replaced  by  assigning  the  constant 
amplitude  previously  found  for  all  samples.)  The  integration  of  the 
frequencies  to  obtain  phase  is  accomplished  by  simply  adding  the 
new  sample  frequency  to  a  total  every  time  through  the  loop.  This 
method  is  a  reasonably  accurate  numerical  integration  technique  since  the 
doppler  frequency  curve  is  highly  over-sampled  with  respect  to  its  fre¬ 
quency  content.  If  the  user  has  specified  that  the  previously  generated 
reference  function  is  to  be  used,  (i.e.,  LFOCUS=3),  the  portion  of  the 
do-loop  used  for  evaluating  the  instantaneous  reference  phase  and 
the  resultant  complex  reference  function  value  is  bypassed,  thereby 
saving  computation  time  and  allowing  the  resident  reference  values  to 
remain  unaltered. 

An  additional  do-loop  is  required  to  offset  the  signal  phase  such  that 
the  phase  at  midbeam  is  zero.  Prior  to  the  integration  procedure  one  has 
no  knowledge  of  the  value  of  the  integral  at  midbeam  and  thus  must  await 
the  completion  of  the  integration  process  to  determine  the  required  phase 
offset.  The  additional  do-loop  is  used  to  subtract  the  required  phase 
offset  and  to  generate  the  resultant  real  signal  values  (zone  plate  samples) . 

The  next  section  of  the  program  prints  the  Chebyshev  coefficients 
and  lists  the  original  doppler  history  samples  if  the  print  option  has 
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been  specified.  A  statement  causing  ail  interpolated  quantities  to  be 
printed  has  been  overridden  at  the  time  of  this  writing,  due  to  the 
huge  amount  of  data  generated  in  high  time-bandwidth  product  cir¬ 
cumstances.  Following  the  print  section  is  a  portion  of  the  program 
responsible  for  generation  of  the  printer  and  Calcomp  plots  of  signal 
amplitude,  doppler  frequency,  signal  phase,  zone  plate  section,  and 
reference  phase  which  have  been  requested  by  the  user  as  I/O  options. 

The  program  is  now  prepared  to  perform  the  convolution  required 
for  target  signature  compression.  To  insure  an  aperiodic  convolution  zeros 
are  padded  onto  the  arrays  containing  the  signal  and  reference  function. 

(If  LFOCUS=3  is  specified,  all  references  to  the  reference  function  are 
ignored.)  Subroutine  RHARM  is  used  to  Fourier  transform  the  real 
valued  signal.  Subroutine  HARM  is  used  to  transform  the  complex 
reference  function  array.  The  resultant  two  transforms  are  multiplied 
point  by  point,  except  that  all  negative  frequencies  are  set  to  zero. 

(Setting  che  negative  frequencies  to  zero  insures  complex  processing  and 
prohibits  aliasing  due  to  conjugate  images  being  formed.)  Taking  the 
inverse  transform  results  in  the  amplitude  response  of  the  simulated 
SAR  to  the  point  target  excitation.  The  modulus  of  each  point  in  this 
array  is  taken  to  produce  a  real  response.  Following  this  step  the 
response  is  plotted  either  directly  as  an  amplitude  response,  or  as 
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a  log-intensity  plot,  and/or  as  an  expanded  iog  intensity  plot,  as 


specified  by  the  user. 


The  section  of  the  subroutine  PROCES  which  will  generate  the  MTF 
and  return  it  in  a  calling-program  supplied  array  will  follow  next  in  ; 


the  final  version  of  PROCES.  The  exact  interface  requirements  for  the 


image  quality  assessment  routines  are  not  complete  at  this  writing  and 


thus  the  MTF  code  is  not  final.  Basically  the  process  will  be  to  square 


the  amplitude  response  and  then  take  an  additional  FFT.  The  MTF  is 


the  modulus  of  the  resulting  transform. 


The  final  portion  of  subroutine  PROCES  is  dedicated  to  generating 


and  plotting  of  the  ideal  system  response  (i.e. ,  perfectly  matched  filter). 


:  If  the  user  requests  this  option,  the  program  proceeds  to  take  the  modulus  > 

\  I 

i  of  the  reference  function  transform,  filter  out  all  negative  frequencies  to  j 


remove  the  conjugate  image  and  then  square  all  terms.  An  inverse  trans¬ 


form  (via  HARM)  results  in  the  equivalent  output  of  a  perfectly  matched 


filter.  Plots  are  made  of  this  ideal  response  according  to  user  specifi¬ 


cations.  After  this  is  complete  control  is  returned  to  the  calling  program. 


7.4  Processing  Examples 


Examples  of  point  target  processing  will  be  given  in  this  section 


to  illustrate  the  plot  options  currently  offered  and  to  demonstrate  the 


ability  of  the  overall  model  to  determine  the  response  of  an  airborne 
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SAR  sys'o.n  in  the  presence  of  realistic  errors  and/or  radar  parameter 
variations.  Preliminary  plotting  examples  of  each  signal  of  interest  to 
the  user  will  be  shown.  Then,  the  results  of  simulated  system  flights 
will  be  reviewed  which  show  the  effects  of  motion-errors  and  antenna 
amplitude  weighting  patterns  upon  system  response. 

Plot  Options  -  For  any  given  target  compression,  the  user  has 
the  option  of  obtaining  any  or  all  of  eleven  signal  plots  as  listed  below: 

1.  Signal  amplitude  vs.  spatial  position 

2.  Signal  frequency  vs.  spatial  position 

3.  Signal  phase  vs.  spatial  position 

4.  One-dimensional  zone  plate  reconstruction 

5.  Reference  function  phase 

8.  Point  target  amplitude  response 

7.  Point  target  log-intensity  response  (db) 

8.  Expanded  log-intensity  response  (db) 

9.  Ideal  system  amplitude  response 

10.  Log-intensity  ideal  system  response  (db) 

11.  Expanded  log-intensity  response  (db) 

Figure  23  is  an  example  of  a  signal  amplitude  plot.  The  horizontal  axis, 
which  represents  spatial  position,  is  labeled  in  terms  of  slant  range 
feet.  Since  the  storage  medium  in  the  digital  system  is  an  imaginary 
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entity,  it  can  be  made  as  large  as  desired  -  in  this  case  large  enough  such 
that  the  scale  factor  between  storage  distances  and  slant  range  distances 
is  unity.  The  example  shown  in  Fig.  23  represents  the  case  where  the 
beamwidth  at  the  target  range  was  400  feet  with  unity  gain  over  than  interval. 
The  target  radar  cross  section  also  did  not  change  with  aspect  angle 
such  that  a  constant  amplitude  linear  FM  pulse  is  returned  and  recorded 
in  memory. 

Figure 24  represents  an  interpolated  doppler  history  as  produced 
by  the  model.  Note  that  the  frequency  modulation  is  basically  linear 
except  for  some  perturbations  near  the  end  of  the  history  due  to  target 
motion  while  in  the  beam.  Note  that  in  this  example  the  beamwidth  at 
the  target  range  was  200  feet.  The  instantaneous  signal  phase  resulting 
from  integrating  the  frequency  history  of  Fig.  24  is  shown  in  Fig.  25.  . 

Note  that  the  phase  has  been  automatically  offset  such  that  it  is  identically 
zero  at  the  midbeam  position.  The  frequency  perturbations  resulting 
from  target  motion  cause  the  instantaneous  phase  to  be  only  approximately 
quadratic.  The  signal  pulse  is  formed  by  merely  taking  the  cosine  of  the 
instantaneous  signal  phase.  The  resultant  pulse  may  be  thought  of  as  a 
one-dimensional  zone  plate.  The  pulse  associated  with  the  instantaneous 
phase  of  Fig.  25  appears  in  Fig.  26.  . 


Doppler  Frequency  of  Return  as  a  Function  of  Azimuth 
Position.  Nonlinearities  are  Due  to  Target  Motion. 


uctei 


The  instantaneous  phase  of  the  reference  function  is  analytically 
constructed  as  a  quadratic  as  was  described  in  the  preceding  section 
on  theory  of  operation.  The  user  may  request  that  the  reference  phase 
be  plotted.  An  example  is  the  phase  plot  of  Fig.  27,  which  is  the 
reference  phase  associated  with  the  signal  phase  of  Fig.  25. 

Convolution  of  the  signal  pulse  and  the  reference  function  produces 
the  desired  point  target  response.  Generally  the  resulting  response  has 
a  (sin  x)/x  appearance  but  may  be  non-symmetric  and  altered  appreciably 
depending  upon  the  system  errors  being  modeled.  Amplitude 
log-intensity,  or  expanded  log-intensity  plots  may  be  requested  by  the 
user.  For  example,  the  amplitude  of  the  target  response  associated  with 
the  signal  pulse  of  Fig.  26  is  plotted  in  Fig.  28..  Note  the  non-symmetry 
and  offset  from  zero  caused  by  the  target  motion.  (Perfectly  focused 
targets  should  be  centered  on  the  zero  position.)  For  simulation  of  low 
resolution  targets  (low  time-bandwidth  product)  the  sidelobes  of  the 
response  can  be  clearly  seen  in  the  normal  amplitude  or  log-intensity 
plots  (log-intensity  =  10  log10  of  amplitude  squared)  but  for  high  resolution 
targets  the  plots  must  be  expanded  in  order  to  see  the  sidelobe  shapes 
and  levels  clearly.  As  an  example,  consider  the  three  plots  appearing 
in  Fig.  29.  This  illustration  shows  the  three  plot  options  as  applied 
to  a  high  resolution  simulation.  The  top  plot  is  an  amplitude  display 
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Figure  28.  Point  Target  Response  (0.0  Azimuth  Position  Represents 


of  point  target  response.  The  middle  plot  is  simply  10  log1Q  of 
amplitude  squared  with  all  portions  less  than  -40  db  set  to  -40  db. 

Note  that  even  in  the  log  plot  the  sidelobes  are  still  not  clearly  dis¬ 
tinguishable.  The  bottom  plot  represents  an  expanded  log  plot  (labeled 
"Truncated  Log  Response")  in  which  those  portions  of  the  response 
greater  than  -40  db  are  expanded  over  the  full  plot  width.  The  side- 
lobe  levels  and  shapes  become  clearly  discernable  in  the  expanded 
plot.  A  note  of  caution  is  in  order  concerning  the  expanded  log  plots. 

Since  only  the  portion  of  the  response  above  -40  db  is  expanded  over  the 
plot  range,  the  expansion  factor  (hence  the  x-axis  scale  factor)  is  thus  a 
function  of  the  response  itself.  Therefore,  in  comparing  the  expanded 
log  plots  of  two  or  more  point  targets,  one  must  pay  special  attention 
to  the  x-axis  labeling  such  that  the  varying  scale  factor  is  taken  into 
consideration. 

The  same  three  plot  options  mentioned  above  for  the  actual  target 
response  are  also  available  for  the  idealized  system  response  which  the 
software  is  capable  of  generating  as  a  byproduct  of  the  signal  compression. 
Characteristics  of  the  ideal  response  is  that  it  is  always  centered  at 
and  symmetric  about  the  zero  location.  The  idealized  system  response 
associated  with  the  actual  target  response  of  Fig.  28  is  shown  as  an 
amplitude  response  in  Fig.  30. . 
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7.4.1  Demonstration  of  Motion  Errors.  A  simulated  flight  was 


made  to  demonstrate  that  the  model  could  indeed  determine  the  effects 
of  imperfect  motion  compensation.  Two  passes  were  made  past  an  array 
made  up  of  three  targets;  one  each  at  far  range,  mid-range,  and  near 
range.  On  the  first  pass,  the  aircraft  was  modeled  as  flying  a  perfectly 
straight  line  at  a  constant  velocity  of  600  ft/sec  at  an  altitude  of  approxi¬ 
mately  40,000  feet.  The  second  pass  was  made  identical  in  every  respect, 

.  2 

except  that  a  constant  vertical  acceleration  of  3.22  ft/sec  (0.1  g)  was 
applied  to  the  aircraft.  During  this  second  pass  the  simulated  motion 
compensation  system  was  instructed  to  correct  for  motion  using  the  mid¬ 
range  point  as  a  reference. 

The  beamwidth  specified  for  this  example  was  C.02  radians.  The 
mid-range  target  was  at  a  ground  range  of  25  NM  with  the  far  and  near 
range  targets  placed  at  plus  and  minus  5  NM  from  mid-range,  respectively. 
The  antenna  was  squinted  forward  0.01  radians  to  provide  an  azimuth  offset. 
Only  azimuth  processing  was  performed  for  this  example.  Under  these 
conditions,  the  time-bandwidth  product  of  the  resulting  doppler  histories 
(azimuth  TBP)  was  something  over  1,000.  To  insure  no  aliasing  problems 
4096  digital  signal  samples  were  used.  (Whereas,  the  number  of  samples 
required  to  satisfy  the  Nyquist  criterion  is  twice  the  TBP.) 
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The  no-motion  azimuth  responses  for  the  far,  mid,  and  near  range 
targets  are  shown  as  amplitude  and  expanded  log-intensity  plots  in 
Figs.  31,  32,  and  33,  respectively.  In  this  example,  unity  signal 
amplitude  was  assigned  to  all  signal  returns,  thereby  ignoring  the 
effects  of  antenna  pattern  and  propagation  attenuation.  Thus,  the 
magnitude  of  the  far,  mid,  and  near  range  responses  are  all  identical. 

Note  that  for  the  no  motion  case  the  responses  are  for  the  most  part  symmetric, 
and  all  are  situated  at  the  zero  spatial  position,  indicating  that  no  target 
image  position  shift  has  occurred.  For  this  example,  the  autofocus  option 
was  chosen  as  the  focusing  method. 

The  target  responses  for  the  case  of  a  constant  0. 1  g  vertical  air¬ 
craft  acceleration  are  shown  in  Figs.  34,  35,  and  36.  Note  that 
there  is  no  significant  change  in  the  shape  of  target  response  for  mid¬ 
range  (Fig.  35)  from  that  demonstrated  for  no  motion  error.  This 
phenomena  is  as  expected  since  the  on-board  motion  compensation  system 
was  instructed  to  correct  for  motion  with  respect  to  the  mid-range  position. 

The  compensation  method  specified  was  to  correct  for  up-beam  velocity  and 
the  compensation  was  to  be  ideal.  That  is,  the  aircraft  navigation  system 
was  given  the  ability  to  exactly  measure  extraneous  aircraft  accelerations. 

Had  the  system  been  more  realistically  specified  then  the  mid-range 
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Figure  32.  Mid  Range,  No  Motion  Point  Target  Response 
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Figure  34.  Far  Range  Point  Target  Response  with  Motion  Error. 
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Figure  35.  Mid  Range  Point  Target  Response  with  Motion  Error. 
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Figure  36.  Near  Range  Point  Target  Response  with  Motion  Error. 
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target  response  may  have  shown  some  smearing  effects.  The  translation 
of  the  mid-range  response  from  the  expected  zero  position  is  due  to  the 
present  implementation  of  electronic  frequency  offset.  The  electronic 
offset  has  been  instructed  to  add  a  frequency  shift  sufficient  to  make 
the  last  in-beam  time  sample  of  doppler  frequency  have  a  value  of 
0.05  cycles  per  foot  or  above.  In  the  present  case,  this  implied  that 
the  resultant  doppler  histories  had  the  correct  FM  rate,  but  an  additional 
frequency  offset.  Thus,  the  azimuth  position  of  the  target  response 
undergoes  an  apparent  shift. 

The  joint  effects  of  aircraft  motion  and  system  attempts  at  motion 
compensation  are  evident  in  Figs.  34  and  36  corresponding  to 
far  and  near  range  targets,  respectively.  Due  to  extraneous  aircraft 
motion  the  linear  FM  rate  and  mid -beam  frequencies  of  the  targets  are 
no  longer  correctly  matched  to  the  reference  function  associated  with 
each  range.  (To  process  the  targets  under  the  influence  of  motion,  the 

t 

reference  functions  derived  in  the  no-motion  case  were  retained  and 
used  again  for  the  associated  targets  in  the  motion  error  examples.) 

Thus,  the  target  responses  at  far  and  near  range  become  smeared, 
asymmetrical,  and  displaced  from  the  expected  location.  The  motion 
compensation  system  modeled  in  this  run  is  such  that  perfect  compensation 
is  possible  only  at  a  specified  range.  Elsewhere  in  the  range  swath  non¬ 
perfect  compensation  is  the  result. 
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The  purpose  of  this  motion  compensation  example  was  to  provide 
test  runs  of  the  system  software  and  to  demonstrate  that  aircraft  motion 
within  our  model  does  manifest  itself  in  the  resulting  target  responses. 

No  attempt  was  made  to  examine  in  detail  the  responses  under  these 
conditions  and  the  effect  of  such  responses  upon  image  quality.  Such 
an  examination  is  not  within  the  scope  of  the  current  program. 

7.4.2  Demonstration  of  Amplitude  Weighting  Effects  -  Antenna. 

Three  additional  simulated  flights  were  made  to  test  the  system’s  ability 
to  demonstrate  the  effects  of  antenna  pattern  and  amplitude  weighting 
upon  the  azimuth  response  of  point  targets. 

Two  flights  were  made  past  an  array  of  three  targets  situated  at 
far,  mid,  and  near-range.  The  two  flights  were  identical  in  every 
respect  except  that  antenna  amplitude  weighting  was  applied  on  the  second 
pass.  In  this  example,  the  vehicle  was  flying  300  ft/sec  at  an  altitude 
of  5,000  ft  and  the  beamwidth  was  0.02  radians.  The  far,  mid,  and 
near-range  targets  were  situated  at  slant  range  distances  of  25K,  23K, 
and  2 IK  feet,  respectively.  The  swath  width  was  approximately  6K  feet. 
Figures  37  and  38  illustrate  the  amplitude  weighting  associated  with 
pass  one  and  two,  respectively.  The  weighting  displayed  in  Fig.  38  is  of 
the  (sin  x)/x  form  out  to  the  -3  db  point  (approximately  -  —  <x<  ^). 

u  u 
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Uniform  Amplitude  Weigl 


The  target  responses  associated  with  the  uniform  weighting 


situation  are  shown  in  Fig.  39,  40,  and 41  for  the  far,  mid,  and 
near-range  targets  respectively.  The  important  items  to  notice  are 
the  sidelobe  levels  and  relative  response  of  the  three  targets.  As 
expected  for  the  uniform  weighting  case,  the  sidelobes  are  down 
approximately  13  db  for  all  three  targets.  In  this  examole,  the 
theoretically  "correct"  reference  function  was  used.  A  slight  mis-match 
between  signal  and  ideal  reference  due  to  range  walk  and  other  higher  order 
effects  causes  a  slight  amount  of  asymmetry  in  the  target  response. 

(The  gross  asymmetries  in  the  tails  of  the  response  (below  -30  db)  are 
most  likely  due  to  finite  aperture  effects.)  The  variation  in  relative 
peak  amplitudes  of  the  three  responses  is  due  to  the  amplitude  weighting 
in  the  range  direction  assigned  to  the  antenna.  This  range  direction 
amplitude  weighting  is  specified  by  the  user  via  a  table  of  values  in 
IQMOD  (the  simulation  section).  Thus  the  antenna  pattern  for  the  first 
flight  had  uniform  weighting  in  the  azimuth  direction  across  the  beam- 
width  but  a  nonuniform  pattern  in  the  range  direction  and  hence  the 
variation  in  the  amplitude  of  returns. 

The  target  responses  in  Figs.  42,  43,  and  44  are  the  result  of 
applying  the  weighting  shown  in  Fig.  39  to  the  far,  mid,  and  near¬ 
range  target  returns.  As  before  there  is  a  relative  amplitude  variance 
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Figure  40*  Uniform  Amplitude  Weighting  Mid  Range  Response. 
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Figure  41.  Uniform  Amplitude  Weighting  Near  Range  Response 
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Figure  42.  Non-Uniform  Weighting  Far  Range  Response 
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(a)  Amplitude  Response 


Figure  43.  Non-Uniform  Weighting  Mid  Range  Response 
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between  the  three  responses  which  is  due  to  the  antenna  pattern  in  the 
range  direction.  The  effect  of  the  amplitude  weighting  in  the  azimuth 
direction  is  shown  in  the  response  sidelobes.  They  are  now  approximately 
23  db  down  as  compared  to  13  db  down  with  unity  weighting.  As  expected 
the  point  target  responses  are  basically  the  same  as  for  the  uniform 
amplitude  case  except  that  the  sidelobe  levels  have  been  reduced.  In 
addition,  one  notices  a  slight  broadening  of  the  mainlobe  as  a  result  of 
the  amplitude  weighting. 

One  additional  flight  was  vnade  to  demonstrate  the  flexibility  with 
which  the  antenna  weighting  can  be  applied.  In  this  final  example,  the 
antenna  pattern  was  broadened  to  include  three  sidelobes  of  the  (sin  x)/x 
pattern  as  shown  in  the  linear  FM  puloe  of  Fig.  45.  To  be  able  to 
process  the  required  time-bandwidth  product  associated  with  a  full 
antenna  pattern  signal,  the  -3  db  beamwidth  had  to  be  halved  from  the 
previous  example.  Therefore,  the  responses  associated  with  the  signal 
pulse  of  Fig.  45  should  not  be  compared  with  any  previous  responses. 
This  example  is  included  solely  to  demonstrate  the  flexibility  of  generating 
arbitrary  antenna  patterns.  The  far,  mid,  and  near-range  responses 
are  shown  in  Figs.  46,  47  and  48  respectively.  The  flat- top  on 
the  mainlobe  of  each  target  response  is  indicative  of  the  mainlobe 
broadening  which  generally  occurs  with  antenna  weighting  schemes. 
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Figure  45.  Three  Sidelobe  Sin  x/x  Antenna  Weighted  Doppler  History. 
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Figure  46.  Far  Range  Target  Response 
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Figure  47.  Mid  Range  Target  Response 
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(a)  Amplitude  Response 


(b)  Expanded  Log-Intensity  Response 
Figure  48»  Near  Range  Target  Response 
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The  sidelobes  of  the  responses  associated  with  this  antenna  pattern 
are  down  about  17  db. 


7 . 5  Summan 


The  impulse  response  capability  has  bcou  included  in  the  model 
to  provide  quantitative  imput  to  the  image  quality  analysis  effort  to  be 
described  in  the  next  section.  The  previous  examples  have  been  included 
to  demonstrate  the  flexibility  and  power  of  the  technique.  It  should  be 
noted  that  we  have  constructed  3.  one  dimensional  routine  equally  applicable 
in  range  or  azimuth  for  the  explicit  purpose  of  providing  detailed  quanti¬ 
tative  data  for  SAR  system  image  quality  analysis. 
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APPENDIX 


Flow  Diagrams  for  the  Impulse  Response  Software 

This  appendix  contains  flow  diagrams  and  listings  for  the  software 
module  described  in  the  section  of  this  report  entitled,  "Impulse 
Response  Generation. "  The  theory  of  operation,  a  narrative  description, 
and  some  processing  examples  are  all  contained  within  that  section. 
Impulse  Response  Flow  Diagram 

The  symbolism  in  the  following  diagrams  is  defined  as  follows: 

-  an  off-page  connector 


-  a  decision  or  branch  point 


-  processing  or  operation  block 


-  a  subroutine  call  (usually  accompanied 
by  an  explanatory  note) 


5' 


-  an  input  or  output  operation 


/ 

L 


-  represents  a  do-loop 


-  entry  point  and  return  points 


t 


-  assignment  indicator 

t 


Software  Module  Listings 

The  listings  of  the  Fortran  source  code  associated  with  the  impulse 
response  module  (IRMOD)  follow  the  flow  diagram. 
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Note  1:  I/O  control  parameter  controls 
the  printing  and  plotting  of 
all  optional  I/O  data 


Note  2:  Chebyshev  curve  fit  to 

signal  Joppler  frequencies. 
Coefficients  used  later  to 
obtain  equidistant  film- 
space  samples. 
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Note  3:  Chebyshev  curve  fit  to  signal 

Amplitudes.  Coefficients  used  later 
to  obtain  equidistant  fils-space 
samples. 


UNIFORM  YES 
AMPLITUDE/ 


CALL  SFIT 


(Note  3) 


IER  *  2 


JT 


AMPCHK 


AMPCHK 


.FALSE. 


.TRUE. 


ALPHA  «■ 
C0EFF(2)/INT 
(Note  4) 


Note  4:  ALPHA  equals  one-half  the  linear  FM 
rate.  Here  ALPHA  is  defaulting  to 
Chebyshev  least  squares  linear  FM 
rate  using  previous  curve  fit  coefficients 


LFOCUS  «  1 


ALPHA  * 
FMRATE/2 
(Note  5) 


Note  5:  For  theoretical  focus  option  ALPHA 
is  reassigned  the  user-supplier 
linear  FM  rate. 


F 


Note  6: 


The  mid  beam  doppler  frequency 
defaults  to  that  actually  observed 
during  target  flyby. 


For  theoretical  focus  option  the 
midbeam  frequency  is  reassigned  to 
the  user  supplied  value. 


CNPS  uses  sample  position  and 
Chebyshev  coefficients  to  deter¬ 
mine  instantaneous  doppler 
frequency  at  that  point. 


YES 


/ PRINT  THE 
I  FREQUENCY 
/  CHEBYSHEV 

'coefficients 


YES 


PRINT  THE 
/AMPLITUDE 
'CHEBYSHEV 
/COEFFICIENT* 


PRINT 

UNIFORM 

AMPLITUDE 

MESSAGE 


PRINT  ORIG/ 
POSITIONS 
' AMPLITUDES j 
/&  PHASES 


I 


PRINT 
INTER¬ 
POLATED 
DATA 


-163- 


-164- 


N 


YES 


-166- 


Note  12: 


•] 
v 


▼ . _ 

■  j  i 

it  CALL  RHARM 

;i 

i|  (Note  12) 

.  ! 


I 

,  -  _ 

j  MULTIPLY 
j  TRANSFORMS 


TOGETHER 
(Note  14) 

ZERO  ALL 
NEGATIVE 


FREQUENCIES 
^  (Note  15) 


lj  ! 

j'  CALL  H^RM  | 
I;  (Note  16)  j 


Note  13: 


Note  14 : 


Note  15: 


Note  16: 


Real  Fourier  transform  of 
zone  plate  section  (one¬ 
dimensional)  in  preparation 
for  convolution. 


Complex  Fourier  transform 
of  reference  function  in 
preparation  for  convolution. 


Multiply  positive  frequencies 
of  transforms  together  to 
achieve  complex  convolution 
of  reference  and  radar  signals. 

Those  negative  frequencies 
resulting  from  real  transform 
of  signal  are  set  to  zero  to 
Insure  complex  processing. 


Inverse  transform  to  complete 
the  convolution. 


TAKE 

MAGNITUDE 

OF 

RESPONSE 


PRINTER 
PLOT  OF 
TARGET 
RESPONSE  / 


CALCOMP 
PLOT  OF 
TARGET 
RESPONSE 


t 


CALCULATE 

Ml'F 

(Note  19) 


YES 


Note  19; 


CALL 
NRMLOG 
(Note  17) 


Note  17 : 


r 


▼ 


The  Modulation  Transfer  Function 
is  the  Modulus  of  the  Fourier 
transform  of  the  intensity 
impulse  response.  The  MTF  is 
generated  and  returned  in  an 
array  to  the  main  program. 

NRMLOG  clips  the  response 
below  a  specified  threshold, 
takes  20  times  the  log  of 
each  point,  and  finds  the 
locations  of  maximum  value 
and  1st  point  above  threshold. 
The  data  is  normalized  such 
that  maximum  is  0.0  db. 
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Note  18: 


To  get  ideal  response  one 
convolves  reference  function 
with  itself.  This  is  equivalent 
to  taking  the  modulus  squared 
of  the  reference  function 
transform  and  performing  an 
inverse  transform. 
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